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Abstract Accreting Millisecond X-Ray Pulsars (AMXPs) are astrophysical labora- 
tories without parallel in the study of extreme physics. In this chapter we review 
the past decade of discoveries in the field. We summarize the observations of the 
fourteen known AMXPs, with a particular emphasis on the multi- wavelength obser- 
vations that have been carried out since the discovery of the first AMXP in 1998. We 
review accretion torque theory, the pulse formation process, and how AMXP obser- 
vations have changed our view on the interaction of plasma and magnetic fields in 
strong gravity. We also explain how the AMXPs have deepened our understanding of 
the thermonuclear burst process, in particular the phenomenon of burst oscillations. 
We conclude with a discussion of the open problems that remain to be addressed in 
the future. 



1 Introduction 

Neutron stars (NSs), amongst the most extreme astrophysical objects in the Uni- 
verse, allow us to study physics in regimes that cannot be accessed by terrestrial 
laboratories. They play a key role in the study of fundamental problems including 
the equation of state (EoS) of ultra-dense matter, the production of gravitational 
waves, dense matter superfluidity and superconductivity, and the generation and 
evolution of ultra-strong magnetic fields. Since the discovery of NSs as radio pulsars 
in 1967, many different classes have been discovered including more than ~ 130 
NSs in low mass X-ray binaries (LMXBs). In LMXBs the NS accretes matter from 
a non-collapsed stellar companion (with mass M < 1 M Q ) via an accretion disk. This 
chapter focuses on a subgroup of the LMXBs, the accreting millisecond X-ray pul- 
sars (AMXPs). In the AMXPs, the gas stripped from the companion is channeled 
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out of the accretion disk and onto the magnetic poles of the rotating neutron star, 
giving rise to X-ray pulsations at the spin frequency. We will explore the details of 
this process, why it is so rare (the AMXPs are a small class), the physics of the disk- 
magnetosphere interaction, and how AMXPs can be used to probe extreme physics. 

Immediately after the discovery of the first millisecond radio pulsar in 1982 (TJ, 
LMXBs were identified as possible incubators for millisecond pulsars. It was sug- 
gested that LMXBs might be responsible for the conversion of slow NSs with high 
magnetic field (B ~ 10 12 G), into a rapidly spinning objects with a relatively weak 
magnetic field (B ~ 10 8 G). Two independent papers published in 1982 [2 , 3 1 pro- 
posed transfer of angular momentum through accretion as the mechanism responsi- 
ble for the spin-up of pulsars. This is known as the recycling scenario (see [:4| for 
an excellent review). This name originates from the fact that radio pulsars switch 
off their pulsed radio emission after entering the so-called "pulsar graveyard". If 
this happens while the NS is in a binary with a non-collapsed low or intermedi- 
ate mass stellar companion, binary evolution of the system [5 , 6] can bring the 
companion into Roche lobe contact and trigger a prolonged epoch of mass trans- 
fer from the companion (donor) towards the NS (accretor). The mass is transferred 
with large specific angular momentum and the NS is spun-up by the resulting ac- 
cretion torques. Once the mass transfer episode terminates, the NS might eventually 
switch on again as a "recycled" millisecond radio pulsar. 

The first AMXP (SAX J1808.4-3658), found in 1998 with the Rossi X-ray Tim- 
ing Explorer (RXTE) satellite [7| provided a beautiful confirmation of the recycling 
scenario. Thirteen more AMXPs have since been found, with spin frequencies from 
182 to 599 Hz. Another important milestone came with the discovery of a radio mil- 
lisecond pulsar in a binary system for which archival optical observations, taken ~ 
10 years before the radio pulsar discovery, show evidence for an accretion disk [8 1. 
This may be the first NS observed to have switched on as a radio pulsar after being 
an X-ray binary. 

The RXTE observatory has played an extraordinary role by discovering many 
systems of this kind and by collecting extensive data records of each outburst de- 
tected during its fifteen year lifetime. The excellent timing capabilities of RXTE 
have brought new means to study NSs with coherent X-ray timing, and helped to 
constrain the long term properties of many AMXPs over a baseline of more than a 
decade. Observation of the orbital Doppler shift of the AMXP pulse frequency con- 
tains information on the orbital parameters of the binary and their evolution in time. 
Binary evolution has benefited from the study of AMXPs [9 10 11] which are now 
known to include ultra-compact systems (orbital period Pf, < 80 min) with white 
dwarf companions, compact systems (Pb — 1.5 — 3 hr) with brown dwarf donors and 
wider systems (Pb ~ 3 .5 — 20 hr) with main sequence stars. Other X-ray and gamma 
ray space missions like XMM-Newton, INTEGRAL, Chandra, Swift and HETE have 
also played an important role in discovering and understanding the spectral and 
timing properties of these objects. Multiwavelength observations covering radio, in- 
frared, optical and UV wavelengths have also illuminated different aspects of these 
fascinating systems. Several optical and infrared counterparts have been identified 
with ground based observations and in some cases have led to the discovery of the 
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spectral type of the donor, while radio and infrared observations have revealed the 
possible presence of jets. 

This chapter is structured as follows: 

Section 2. An overview of the AMXP family. 
Section 3. Observations of individual AMXPs 
Section 4. Coherent timing analysis and accretion torques. 
Section 5. X-ray pulse properties, their formation and use. 
Section 6. Long-term evolution of spins and orbital parameters. 
Section 7. Thermonuclear bursts and burst oscillations. 
Section 8. Aperiodic phenomena including kilohertz QPOs. 
Section 9. Future developments and open questions. 



2 The Accreting Millisecond X-ray Pulsar Family 

An AMXP is an accretion powered X-ray pulsar spinning at frequencies v > 100 Hz 
bound in a binary system with a donor companion of mass M < IM®. All known 
AMXPs have donor stars that transfer mass via Roche lobe overflow (RLOF). This 
definition excludes binary millisecond pulsars, which are detached systems with a 
pulsar powered by its rotational energy, X-ray pulsars in high mass X-ray binaries, 
symbiotic X-ray binaries and the slow X-ray pulsars in LMXBs with v < 100 Hz. 
Several LMXBs show millisecond oscillations during thermonuclear bursts that ig- 
nite on their surface. These nuclear powered X-ray pulsars (NXPs), or burst oscilla- 
tion sources, are not AMXPs, since they are powered by nuclear burning rather than 
channeled accretion (although note that some AMXPs also show burst oscillations). 

Table Q] reports the main characteristics of the 22 known pulsating NSs in 
LMXBs: 14 AMXPs, 3 slow accreting X-ray pulsars in LMXBs, 1 slow accret- 
ing X-ray pulsar in an intermediate mass X-ray binary (IMXB), 3 symbiotic X-ray 
binaries accreting wind from a K/M giant companion (compatible with being of low 
mass type) and 1 mildly recycled accreting pulsar recently discovered in the glob- 
ular cluster Terzan 5. This latter source is an 11 Hz X-ray pulsar with a magnetic 
field of 10 9 — 10 10 G and is of particular interest because it might be the only known 
accreting pulsar in the process of becoming an AMXP on a short timescale lfT2l . 
The other four slow pulsars are very different systems, in the sense that they have 
most likely followed a completely different evolutionary history, have strong mag- 
netic fields (B ~ 10 12 G) and may never reach millisecond periods. Symbiotic X-ray 
binaries are also different from the AMXPs since they are wide binaries (Pf, > 50 
days) and their NSs have extremely long spin periods caused by a prolonged phase 
of spin-down during a wind accretion process. 

All of the AMXPs are transient systems, with accretion disks that run through 
cycles of outburst and quiescence. X-ray pulsations have never been observed in 
quiescence, although the upper limits are unconstraining due to the poor photon 
flux. The shortest outburst recurrence time is one month for the globular cluster 
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source NGC 6440 X-2, with an outburst duration of less than 4-5 days, whereas the 
longest outburst, from HETE J1900. 1-2455, has so far lasted for almost 7 years. 
The small group of AMXPs therefore form a rather heterogeneous class of objects. 
There are, however, some common features shared by all AMXPs: 



• Outburst luminosities are usually faint, suggesting low time averaged mass 
accretion rates (Section|3]l. 

• Spin frequencies appear to be uniformly distributed with an abrupt cutoff 
at about ~ 700 Hz (Section©. 

• Ultra-compact binaries are rather common, comprising about 40% of the 
total AMXP population. 

• Very small donors are preferred, with masses almost always below O.2M . 



Table 1 Accreting X-ray Pulsars in Low Mass X-ray Binaries 



Source 


v s 

(Hz) 


Pi 
(hr) 


fx 

(**o) 




M c „„■„ 

w4) 


Companion Type Ref. 




Accreting Millisecond Pulsars 


SAX J1808.4-3658 


401 


2.01 


3.8 x 10- 


-j 


0.043 


BD 




14 


1 1511 1611 171 


XTE J1751-305 


435 


0.71 


1.3 x 10 


-6 


0.014 


He WD 


m 


19 




XTE J0929-314 


185 


0.73 


2.9 x 10 


-7 


0.0083 


He WD 


[20H21 




XTE J 807-294 


190 


0.67 


1.5 x 10 


-7 


0.0066 


C/O WD 


(22 


23 




XTE J1814-338 


314 


4.27 


2.0 x 10 


-3 


0.17 


MS 


[27 


28 




IGR J0029 1+5934 


599 


2.46 


2.8 x 10 


-5 


0.039 


BD 


(29 


30 




HETE J1900.1-2455 


377 


1.39 


2.0 x 10- 


-6 


0.016 


BD 


[35] 




SwiftJ 1756.9-2508 


182 


0.91 


1.6x 10" 


-7 


0.007 


He WD 


(36 


37 




Aql X-l 


550 


18.95 


N/A 




0.6° 


MS 


38 


39 






















SAXJ1748.9-2021 


442 


8.77 


4.8 x 10 


-4 


0.1 


MS/SubG ? 


|40j 


141] 


NGC6440 X-2 


206 


0.95 


1.6x 10- 


-7 


0.0067 


He WD 


m 




IGR J 175 11-3057 


245 


3.47 


1.1 x 10 


-3 


0.13 


MS 


[43] 


|44] 


Swift J 1749.4-2807 


518 


8.82 


5.5 x 10- 


2 


0.59 


MS 


EH 




IGR J 17498-2921 


401 


3.84 


2.0 x 10 


-3 


0.17 


MS 


|46| 




Mildly Recycled X-ray Pulsar 


IGR 117480-2466 


11 


21.27 


2.1 x 10 


-2 


0.4 


SubG 


47 48. 49 


Slow Pulsars in LMXBs & IMXBs 


2A 1822-371 


1.7 


5.57 


2x 10~ 2 




0.33 




[50] 




4U 1626-67 


0.13 


0.69 


1.3 x 10 


-6 


0.06 


WD or He star 


I51II52II53I 


GRO 1744-28 


2.14 


282.24 


1.3 x 10 


-4 


<0.3 


Giant 


[5l]|54][55] 


Her X-l 


0.81 


40.80 


0.85 




~2 D 


MS 


|51 


53 




Symbiotic X-ray Binaries 


GX 1+4 


6.3 x 10 


1 1161 days 








M5 Giant 


[54 


56 


1571 


4U 1954+31 


5.5 x 10 


-5 








M4 Giant 


m 


57 




IGR J16358^t724 4 


1.7 x 10 


-4 








K/M Giant 


58 57 





Vj is the spin frequency, P% the orbital period, f x is the X-ray mass function, M c mi „ is the minimum companion mass 
for an assumed NS mass of IAMq. The companion types are: WD = White Dwarf, BD= Brown Dwarf, MS = Main 
Sequence, SubG = Sub-Giant, He Core = Helium Star. 

" The donor mass is inferred from photometric data and does represent the most likely mass. 
b Binaiy with parameters that are still compatible with an intermediate/high mass donor. 
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2.1 Intermittency 

Until 2007 it was believed that AMXPs showed X-ray pulsations throughout out- 
bursts. The seventh AMXP (HETE J 1900. 1-2455), however, showed an unexpected 
new behaviour [61|. For the first w 20 days of its outburst (which started in 2005) 
it showed typical AMXP pulsations. The pulsations then became intermittent, ap- 
pearing and disappearing on different timescales for the next 2.5 years. Pulsation 
disappeared and never reappeared again after MJD 54,499 [62 1, with the most strin- 
gent upper limits of 0.05% rms on the fractional pulse amplitude (i.e., a sinusoidal 
fractional amplitude of 0.07%, see Eq.[T4"l>. 

This discovery was exciting, since it may help to bridge the gap between non- 
pulsating LMXBs and AMXPs. It also became immediately clear that other LMXBs 
might show sporadic episodes of pulsations during their outbursts. This has now 
been found to be the case for two other sources: Aql X-l and SAX J1748. 9-2021. 
Aql X-l is perhaps the most striking case recorded so far: coherent pulsations were 
discovered in 1998 RXTE archival data and appeared in only one 120s data segment 
out of a total exposure time of 1.5 Ms from more than 10 years of observations. 
The extremely short pulse episode has raised discussions about whether it really 
originated from magnetic channeled accretion, and whether Aql X-l can truly be 
considered an AMXP. However, the high coherence of the signal leaves little doubt 
about the presence of pulsations, and the accretion-powered origin appears the most 
promising explanation |38l . The case of SAX J 1748. 9-2021 is slightly different: 
pulsations were detected sporadically in several data segments and in three (2001, 
2005 and 2009-2010 EU gOJ ED) out of four outbursts observed (the first being in 
1998). It is unclear why these three systems show pulsations intermittently. In HETE 
J 1900. 1-2455, an increase in pulse fractional amplitude was reported approximately 
in coincidence with the occurrence of Type I X-ray bursts ll6Tl . followed by a steady 
decrease. On other occasions the pulsations appeared a few hours before or after a 
burst, indicating that pulsations might be linked somehow with some yet to be iden- 
tified property of the NS envelope. In SAX J1748. 9-2021 the pulsed amplitudes 
showed some abrupt changes in amplitude and/or phase in coincidence with about 
30% of the observed Type I X-ray bursts [HTI |401 . The pulsations, however, dis- 
played a more diverse behaviour than in HETE J1900. 1-2455, without the typical 
steady decrease of fractional amplitudes. A period of global surface activity during 
which both Type I bursts and pulsations are produced might be at the origin of this 
link [40 1 . The single pulsating episode of Aql X-l had instead no clear connection 
with Type I bursts, even though Aql X-l is a bursting LMXB. 

The spin frequency derivative of HETE J 1900. 1-2455 was measured over a base- 
line of 2.5 years, an unprecedented long baseline for an AMXP (whose outbursts 
last usually less than 100 days). This has indirectly provided hints on the physical 
origin of such a period of global surface activity in intermittent sources. The spin 
frequency derivative exhibited an exponential decay in time that was interpreted as 
evidence of the burial of the NS magnetic field ll62l . It was proposed that intermit- 
tency originates because the magnetic field strength drops by almost three orders of 
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magnitude on a timescale of few hundred days (see also Section [631 for a discussion 
of the model). 

One feature that intermittent pulsars share is that the long term average mass 
accretion rate (M) is higher than for the persistent AMXPs and smaller than for 
the bright non-pulsating systems (like Sco X-l and the other Z sources). Calcu- 
lating the precise value of (M) is rather difficult, since it depends on poorly con- 
strained parameters in most LMXBs, like the distance d, the X-ray to bolometric 
flux conversion and the recurrence time of the outburst. However, even considering 
these caveats, it seems clear that at least the brightest systems have never showed 
pulsations. Whether it is the high mass accretion rate that determines the lack of 
pulsations or some other properties shared by the brightest systems is still unclear. 
See Section [6] for an extended discussion on the mechanism that might prevent the 
formation of pulsations in LMXBs. 



3 Observations of the AMXPs 

In this section we discuss the main characteristics of each AMXP, including num- 
ber of outbursts observed, luminosity, distance, orbital parameters and (at the end of 
each subsection) radio/IR/optical counterparts. We devote Section|4]and|6]to the dis- 
cussion of the results on the pulsar rotational parameters and their secular evolution, 
respectively. We briefly mention results relating to X-ray bursts, burst oscillations 
and aperiodic variability but refer to Section|7]and|8]for a general discussion of these 
topics. The astrometric position of each AMXP is given in TablefJ] 

Table 2 Astrometric Position of AMXPs 



Source Right Ascension Declination Error Observatory Reference 





(J2000) 


(J2000) 












SAX J 1808.4-3658 


18:08:27.62 


-36:58:43.3 


0' 


.15 


6.5 m Baade (Magellan I) 


14 




XTE J1751-305 


17:51:13.49(5) 


-30:37:23.4(6) 






Chandra 


18 




XTE J 1807-294 
IGR J00291+5934 


18:06:59.80 
00:29:03.05(1) 


-29:24:30 
+59:34:18.93(5) 




1" 


Chandra 

Multiple Optical Obs. 


22 
64 




XTE J1814-338 


18:13:39.04 


-33:46:22.3 


0> 


.2 (90% 


c.l.) Magellan 


65 




XTE J0929-314 


09:29:20.19 


-31:23:03.2 


0' 


.1 


Mt. Canopus 


66 




Swift J 1756.9-2508 


17:56:57.35 


25:06:27.8 


3' 


.5 


SwiffXRT 


36 




Aql X-l 


19:11:16.0251(3) +00:35:05.891(4) 






EVN 


67 




SAX J1748. 9-2021 


17:48"52.163 


20:21:32.40 


ff 


.6 


Chandra 


68 




NGC6440 X-2 


17:48:52.76(2) 


20:21:24.0(1) 






Chandra 


69 




Swift J 1749.4-2807 


17:49:31.83 


-28:08:04.7 


r 


.6 (90% 


c.l.) Swift/XBJ 


70] 


IGR J 175 11-3057 


17:51:08.66(1) 


-30:57:41.0(1) 


0' 


.6 


Chandra 


71] 


IGR J 17498-2921 


17:49:55.38 


-29:19:19.7 


ry 


.6 


Chandra 


72] 


HETE J 1900. 1-2455 


19:00:08.65 


-24:55:13.7 


0' 


.2 


Palomar 


73] 



The errors in parentheses refer to R.A. and DEC separately, whereas those that appear in the column "Error 1 refer to 
both coordinates. 
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3.1 SAX J 1808.4-3658 

The source SAX J1808. 4-3658 was discovered with the X-ray satellite BeppoSAX 
during an outburst in September 1996 [74|. The source showed via three type I 
X-ray bursts that the compact object in the binary is a NS. The 2-28 keV peak lumi- 
nosity of 3 x 10 36 ergs -1 IT751 is rather faint and below the average peak luminosity 
reached by other LMXBs. X-ray pulsations were not detected during this outburst, 
with poorly constraining upper limits of 20% on the pulsed fraction. Coherent pul- 
sations at 401 Hz were discovered instead in 1998, during the second observed out- 
burst 1 7 1, thanks to the better sensitivity of the Proportional Counter Array (PC A) 
on RXTE. An orbital modulation of 2.01 hr was detected using Doppler delays in the 
coherent timing data [76| and it was suggested that the companion of the pulsar is 
a heated brown dwarf with mass of w 0.05 M Q 1771 . No thermonuclear X-ray bursts 
were observed during the 1998 outburst, but re-analysis of the 1996 data provided 
marginal evidence for burst oscillations at the NS spin frequency ll75ll . 

SAX J 1808.4-3658 went into outburst again in 2000, 2002, 2005, 2008 and 201 1, 
with an approximate recurrence time of about 1.6-3.3 yr, and is the best sampled and 
studied of all AMXPs. All of the outbursts showed faint luminosities, with peaks 
always below 10 37 ergs _1 , even when considering broad energy bands. The 2000 
outburst was poorly sampled due to solar constraints, and was observed only during 
the "flaring tail", a peculiar outburst phase displayed by only a very small number of 
X-ray sources. The typical SAX J1808.4-3658 outburst can be split in five phases: a 
fast rise, with a steep increase in luminosity lasting only a few days, a peak, a slow 
decay stage, a fast decay phase and the flaring tail. The first four phases are typi- 
cal of several X-ray binaries and dwarf novae and can in principle be partially ex- 
plained with the disk instability model. The flaring tail instead shows bumps, called 
"refiares", with a quasi-oscillatory behaviour of a few days and a variation in lumi- 
nosity of up to three orders of magnitude on timescales ~ 1-2 days [78, 79, 80 1 . The 
flaring tail has no clear explanation within the disk instability model [81 1, but has 
been observed in all outbursts [80 82 1. Pulsations are still detected during this low 
luminosity stage, and a strong ~ 1 Hz oscillation is observed to modulate the refiares 
(see Fig. [D [83 80 1). This has been interpreted as a possible signature of the onset 
of a disk instability with radial perturbation of the magnetospheric boundary (the 
so-called "Spruit and Taam instability" [84, 85 1) close to the onset of the"propeller 
stage" 1 80 1 . The 1 Hz modulation, however, has not appeared in all the flaring tails 
of the different outbursts, but only in the years 2000, 2002 and 2005. The presence 
of the 1 Hz modulation during the 1998 flaring tail cannot be excluded, however, 
due to lack of observations. 

The study of the high frequency aperiodic variability of SAX J1808. 4-3658 led 
to a breakthrough when twin kHz Quasi Periodic Oscillations (QPOs) at a frequency 
of approximately 500 and 700 Hz were discovered for the first time in a source 
with a well established spin frequency [86|. The twin kHz QPOs (Section [8) were 
observed when SAX J1808. 4-3658 reached the peak luminosity during the 2002 
outburst, and their frequency separation of 196 ±4 Hz is consistent with being at 
half the spin frequency. This suggested a possible link between spin frequency and 
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Fig. 1 Small portion (15 
s) of the 2005 X-ray light- 
curve of SAX J1808.4-3658, 
showing the 1 Hz flaring 
phenomenon during a reflare. 
The amplitude of the 1 Hz 
oscillation reaches 125% rms 
(Figure from [80 1). 



1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 
Time [s] 



kHz QPOs and hence refuted the beat-frequency model for the formation of kHz 
QPOs. In this model the upper kHz QPO at 700 Hz reflects an orbital frequency 
in the inner accretion disk, whereas the lower kHz QPO should appear as a beat 
frequency between the upper kHz QPO and the spin frequency. In this model the 
twin kHz QPOs must be separated by a frequency equal to the spin frequency and 
not half its value, as observed. The upper kHz QPO was observed during most of the 
outburst, and another QPO at 410 Hz also appeared during a few observations [86|. 
The origin of this latter QPO is unclear and it has been suggested [ 86 1 that it might be 
related to the side-band phenomenon observed in other LMXBs 1 87 1 . In particular, it 
was suggested that the 410 Hz QPO might be a sideband of the pulsation, created by 
a resonance occurring at the radius in the disk where the general relativistic vertical 
epicyclic frequency matches the spin frequency of the pulsar. 

The 3-150 keV X-ray spectrum of the 1998 outburst had a remarkably stable 
power law shape with photon index of ~ 2 and a high energy cutoff at ~ 100 keV 
[88 1. Further spectral analysis provided evidence for a two-component model: a 
blackbody at soft energies and a hard Comptonization component at higher energies 
(89). The blackbody is interpreted as the heated hot spot on the NS surface, whereas 
the Comptonization is produced in the accretion shock created at the bottom of the 
magnetic field lines as the plasma abruptly decelerates close to the NS surface. The 
presence of an accretion disk was detected much later at lower energies, with obser- 
vations taken in 2008 with the EPIC-pn camera on the XMM-Newton telescope. Its 
large sensitivity at soft energies (down to about 0.5 keV) well below the nominal 2 
keV limit of RXTE, allowed the detection of the typical cold accretion disk signa- 
ture at a temperature of 0.2 keV [16. 90 1. The signature of a fluorescent relativistic 
iron Ka emission line profile was also found. A similar result was obtained with 
combined XMM-Newton and Suzaku data ifTTl . Spectral modeling of the iron line 
constrained the magnetic field of the pulsar to be ~ 3 x 10 8 G at the poles lfl7l[T6l . 
Simultaneous spectral modeling of the inner disk radius and pulse profile shapes of 
the 2002 outburst 11911 lead to a similar constraint of the magnetic field (B ~ 10 8 G). 
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a _ b 




10 20 30 40 -0.5 0.0 0.5 1.0 



Seconds since start of burst Seconds since start of burst 

Fig. 2 Left panel: a typical X-ray burst observed in the 2-60 keV energy band on 1 8 October 2002 
with a dynamical power spectrum of barycentered data superimposed. The contours refer to power 
levels of 2 s power spectra spaced by 0.25 s. The burst oscillations are clearly detected during 
the rise, overshooting the spin frequency of SAX J1808. 4-3658 at 401 Hz. The oscillations then 
disappear during the photospheric radius expansion close to the burst peak and reappear in the burst 
tail at a frequency consistent with the 401 Hz of the pulsar. Right panel: a zoom in of the burst 
rise and peak, with finer power spectra resolution of 0.25 s spaced by 0.03125 s. The frequency 
overshoot is clearly visible, with the BO frequency reaching ~ 403 Hz (Figure from 1921 ). 



Thermonuclear bursts were observed in 1996, 2002, 2005, 2008 and 2011 J75J 
|92l [82ll . Most of the bursts exhibit photospheric radius expansion (PRE), where the 
luminosity reaches the Eddington limit, lifting the photosphere off the surface of the 
NS until the flux dies down. Such bursts can be used as standard candles [93|. The 
distance estimated using this method is 2.5-3.6 kpc [82, 75]. Note that a different 
lower limit of 3.4 kpc is reported in [82 1. This lower limit is based on the assumption 
that the long-term mass transfer rate is driven purely by loss of angular momentum 
in the binary via emission of gravitational radiation, which may not be a good ap- 
proximation (see Section^- All bursts observed in the RXTE era have shown burst 
oscillations (with a possible marginal detection in one burst in 1996 observed with 
BeppoSAX) with an amplitude of a few percent rms. SAX J1808. 4-3658 provided 
the first robust confirmation that burst oscillation frequency was, to within a few Hz, 
the spin frequency of the star ([92] and Figure [2). 

An optical/IR counterpart (V4584 Sagittarii) was discovered during the 1998 out- 
burst, coincident with the position of SAX J1808.4-3658 [94|. The reported mag- 
nitudes of the candidate were V = 16.6, R = 16.1, 1 = 15.6, J = 15.0, H = 14.4, K 
= 13.8, with an uncertainty of 0.2 mag in VRI and < 0.1 mag in JHK. The V band 
was further monitored 11951 and a possible sinusoidal modulation at the 2 hr period 
of the binary identified, together with a decay in the luminosity as the outburst pro- 
gressed in its decay stage. A multiband optical/IR photometric study of the optical 
counterpart during the 1998 outburst revealed an optical flux consistent with an X- 
ray heated accretion disk and an inclination of cos; = O.65+0J3 (90%c.l.)|96[.The 
IR observations, however, showed an excess with respect to an accretion disk plus 
irradiated donor star model. This excess was transient in nature, as it was detected 
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only during one observation, the rest being consistent with the irradiated accretion 
disk plus donor star model. Radio observations carried out one week later (April 
27, 1998) with the Australia Telescope Compact Array (ATCA) revealed a transient 
radio counterpart with a flux of ~ 0.8 mJy with no further detection at later epochs 
during the 1998 outburst or during quiescence. This observation was interpreted as 
synchrotron emission which might in principle constitute a viable explanation also 
for the IR transient excess. This suggests that material can be ejected from the binary 
via relativistic jets and/or outflows. 

Optical observations in quiescence have been performed extensively for over a 
decade 1 97, 98 , 99] . In the first observation, a clear optical excess was detected, well 
above the level expected from residual X-ray irradiation of the donor star that is still 
present during quiescence [ 97] . A clear modulation at the orbital period of the binary 
in IR bands was also identified Il98"lll00ll99ll . To explain the observations, a strong 
source of extra irradiation was required. In particular, the initial donor mass estimate 
of 0.05 M [77 1 was revised to 0.07-0.1 1M , to account for the large entropy of the 
donor in quiescence [98 1. The optical excess was first interpreted as evidence for 
the turning on of an active radio pulsar when accretion halts f 1011 11021 . Whether 
this is the correct interpretation is still to be verified, but it appears a plausible ex- 
planation even though no radio pulsations have been detected [ 103 1 and alternative 
possibilities exist (see Section 6 and f\3\ ). 



3.2 XTE J 175 1-305 

The AMXP XTE J175 1-305 was discovered by RXTE on April 3, 2002 and coher- 
ent pulsations at 435 Hz were observed immediately [18]. The orbital period of 42 
minutes makes this AMXP an ultra-compact binary with a heated He or C/O white 
dwarf of minimum donor mass O.O13-O.O17M Q (depending on the NS mass) flTTj . 
Archival observations taken with the All-Sky-Monitor (ASM) aboard RXTE showed 
that a faint outburst had already occurred in 1998. Three other outbursts were ob- 
served in 2005, 2007 and 2009. The first two had very faint peak luminosities of 
about 10%-20% the value reached in 2002. The 2005 outburst in particular was ob- 
served by INTEGRAL [ 105 1 and lasted only 2 days, with only one follow-up RXTE 
observation II 1 061 - No high resolution timing data were taken, and the association of 
the outburst with XTE J1751-305 remains dubious given the large source position 
error of the two X-ray observatories. On April 5, 2007, XTE J 175 1-305 was ob- 
served with RXTE during routine galactic bulge monitoring observations [11071 . No 
pulsations were seen because of the dim flux (10 mCrab in the 2-10 keV) and the 
short bulge scan exposure (~ 50 s). However, a follow-up Swift observation made 
the identification of the source certain thanks to the high angular resolution of the 
XRT telescope [ 108]. The outburst had a fiuence of about 10~ 4 ergcm~ 2 , compara- 
ble to the already low 2002 outburst fiuence of 2.5 x 10~ 3 ergcm~ 2 [18). In 2009 
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INTEGRAL recorded a new outburst 01 091 and RXTE observed it extensively, with 
pulsations immediately detected lll 101 . This outburst has a duration and brightness 
similar to the 2002 event. 

XTE J1751-305 shows no Type I X-ray bursts and its distance is therefore dif- 
ficult to determine. The X-ray spectrum taken with XMM-Newton in 2002 showed 
a power law with spectral index 1 .44 ± 0.01 contributing 83% of the total 0.5-10.0 
keV flux and a blackbody at temperature kT = 1.05 ±0.01 keV 111 II . A broadband 
spectrum analysis of the 2002 outburst was performed by combining RXTE (PCA 
and HEXTE) and XMM-Newton (EPIC-pn) data |[TT2l . The analysis revealed the 
presence of three components, two blackbodies and a hard component, which were 
interpreted as a thermal emission from a cold accretion disk (kT ~ 0.6 keV), a hotter 
(kT ~ 1 keV) hot spot on the NS surface and a thermal Comptonization in plasma 
of temperature kT ~ 40 keV (see Fig |3J. These findings were very similar to the 
spectral modeling described above for SAX J1808. 4-3658. Unsuccessful searches 
for an optical/near IR counterpart were carried out during the 2002 outburst and in 
quiescence 111 131 . However the upper limits are not particularly constraining if one 
assumes that the source is close to the Galactic center (d — 8.5 kpc). 




Fig. 3 Broadband X-ray spectrum of XTE J 175 1-305 obtained by combining XMM-Newton (light 
blue points) and RXTE data (black and dark blue points). Three components are used to model the 
spectrum: a cold blackbody identified with the thermal disk emission (green line), a hot blackbody 
(hot spot; blue line)) and a hard tail created via thermal Comptonization in a plasma (pink line). 
The left spectrum refers to the emission geometry "(a)" in which a small temperature gradient is 
present in the hot spot. The right spectrum uses the "(b)" model, with a large temperature gradient 
between the hot spot and the remaining illuminated NS surface around the accretion shock (image 
taken from 1 1121 ) 
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With a spin of 185 Hz and an orbital period of 43.6 min, XTE J0929-314 was the 
third AMXP to be discovered, and the second in an ultra-compact binary. The source 
was first detected on April 13, 2002 by the RXTE/ASM II 1 1411 . Coherent pulsations 
were first detected in an observation performed on May 2, 2002 [20|. Doppler mod- 
ulation of the pulsations revealed the ultra-compact nature of the binary [21 1. The 
very small mass function (2.7 x 10~ 7 M n ) implied a tiny donor star of minimum 
mass of about 0.01 M Q . This requires a cold He white dwarf companion with a pos- 
sible H-rich envelope that would explain the detection of an H a line at A = 6, 563 A 
111 151 during the 2002 outburst. The 2002 outburst lasted at least 52 days, but only 
34 days were monitored by the high time resolution RXTE/PCA instrument. The 
X-ray spectrum of XTE J0929-314, taken with Chandra, showed a power-law plus 
blackbody similar to XTE J 175 1-305 HE). 

XTE J0929-314 has been detected in several wavebands: radio, near IR and op- 
tical counterparts were identified during the 2002 outburst. An optical counterpart 
consistent with the position of XTE J0929-314 was detected in V, B, R and I bands 
111 171 . Broad band BVRI observations at the Mt Canopus 1-m telescope carried in 
May 2002 revealed a modulation of 10% in the I band with periodicity consistent 
with the orbit of the binary ||66l . An anomalous I-band excess was also observed 
in the first few days of the outburst and was interpreted as synchrotron radiation, 
in addition to the accretion disk/irradiated donor emission. This behavior is similar 
to that reported for SAX J1808.4-3658 in the V and I band. The optical spectrum 
of the source (in outburst) shows a strong emission line at 4,640 A compatible with 
a C/O or a He/N disc HI 181 . Optical spectroscopy also lead to the conclusion that 
the most likely donor of XTE J0929-314 is a C/O white dwarf. A radio counterpart 
(during the outburst) was detected with the Very Large Array (VLA) at 4.86 GHz 
GUI. A source flux of 0.31 ±0.07 mJy (May 3) and 0.36 ±0.05 mJy (May 7) was 
reported from a position only 0.8 arcseconds from the previously detected optical 
position 111 171 . The source has also a very faint optical counterpart (R = 27.2) in 
quiescence [120]. PSF VRI photometry in quiescence [121] also suggested that the 
donor is irradiated by a source which emits in excess (by a factor of « 8) of the 
quiescent X-ray luminosity of XTE J0929-3 14 (L x ~ 4.5 x 10 31 ergs _1 ). 



3.4 XTE J1807-294 

The fourth AMXP was discovered on February 21, 2003 by RXTE 11221 . The source 
is the third ultra-compact AMXP, with the shortest orbital period known to date, 40. 1 
min, and a minimum companion mass of O.OO66M . Simultaneous INTEGRAL, 
XMM-Newton and RXTE spectral analysis [122] constrains the companion mass to 
be less than 0.022M3, implying that the donor is likely a C/O white dwarf rather 
than a He white dwarf, which would require unlikely a priori low inclinations. With 
the exception of the quasi-persistent intermittent AMXP HETE J1900.1-2455, XTE 
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J 1807-294 had the longest AMXP outburst observed so far: w 120 days, of which 
ss 109 days had sufficiently high flux to allow the detection of pulsations. Twin kHz 
QPOs were detected in several observations and the separation A V of the kHz QPOs 
is consistent with the spin frequency of the AMXP [ 123 1 . No counterparts have been 
reported at any wavelength during either outburst or quiescence. A candidate optical 
counterpart detection was reported in 2009 111211 with V ~ 22 . 1 and R ~ 2 1 .4 along 
with a tentative detection in the J band. The counterpart falls within the 1" Chandra 
error circle. No variability has been observed, with upper limits of 0. 1 mag. 



3.5 XTEJ1814-338 

The fifth AMXP was found during routine Galactic center observations by the 
RXTE/PCA on June 5, 2003 |27). The pulse frequency is 314 Hz and the orbital 
period 4.3 hr, giving a minimum companion mass of 0.17M©. The outburst lasted 
for ss 50 days and pulsations were detectable for w 45 days. The source showed 
28 Type I X-ray bursts and RXTE detected burst oscillations at the spin frequency 
in all bursts. The 28th burst showed signs of PRE, allowing a determination of the 
source distance of w 8 kpc. The X-ray position of XTE J18 14-338 is consistent with 
the EXOSAT source EXMS B1810-337, detected in September 1984 E24). If this 
identification is correct, then XTE J18 14-338 has a recurrence time of ~ 20 yr. 

The optical counterpart was first tentatively identified using BVRI data taken 
with the Magellan 6.5 m telescope 11251 . An observation of H and He emission 
lines, with a double-peaked H a line typical of interacting binaries, was reported 
soon after the discovery of the binary [126]. The counterpart was studied in more 
detail and finally unambiguously identified 1631 . The B, V and R counterparts during 
the 2003 outburst had magnitudes between 18 and 19 mag, whereas the I band was 
observed with an excess flux reaching about 17.4 mag. This excess is again similar to 
that reported for SAX J1808.4-3658 and XTE J0929-3 14 and suggests the presence 
of jets and/or outflows. No counterpart was observed in quiescence, with a limiting 
magnitude of R > 23.3 (3(7 c.L, (§3). A faint object with V ~ 23 .3 and R ~ 22.5 was 
identified in ESO Very Large Telescope (VLT) archival data taken in May 20-21, 
2004 11211 . The optical position was found to be consistent with values previously 
reported ll65l and the light curve showed variability consistent with the 4.3 hr orbital 
period of the binary thus confirming that this is indeed the optical counterpart in 
quiescence of XTE J1814-338. 



3.6 IGR J00291+5934 

IGR J0029 1+5934 is the fastest AMXP known, with a spin of 599 Hz. The source 
was discovered by INTEGRAL on 2004 December 2 II 1271 and prolonged RXTE 
observations lead to the discovery of pulsations [ 128 1 and an orbital period of 2.46 
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hr, with a minimum donor mass of about 0.04 M Q 111291 . The 2004 outburst had a 
duration of w 14 days from the INTEGRAL discovery, with a smooth flux decay. On 
2008 August 13 the source started a new outburst [130] and was observed for about 
a week before reaching an X-ray flux level below the sensitivity limit of RXTE. The 
peak flux observed was about half the value reached in 2004. Follow up SwiftfKRT 
observations on August 21 revealed a flux level below 10~ n — 10 _12 ergs _1 and 
an XMM-Newton observation on August 25 provided a 2-10 keV flux of (1.4 ± 
0.3) x 10 _14 ergs _1 II 1 3 1 1 . compatible with the 10~ 13 ergs~' (0.5-10keV) quiescent 
flux level observed two months after the end of the 2004 outburst [ 1 32 1 . On 2008 
September 21, the source re-brightened again, showing another outburst lasting for 
about 14 days. The fluence of the second and third outbursts were very similar and 
it is unclear whether the two 2008 outbursts were distinct [33 1 or part of the same 
outburst. An analysis of archival RXTE/ASM data revealed possible outbursts in 
1998 November and 2001 September, suggesting a recurrence time of 3-4 yr 11331 . 

The pulsar has shown no thermonuclear bursts despite very similar orbital param- 
eters to SAX J1808.4-3658. The donor star of IGR J00291+5934 is also probably 
a X-ray heated brown dwarf as in SAX J1808. 4-3658 ll29l . Broadband spectral ob- 
servations performed during the 2004 outburst with INTEGRAL and RXTE were 
well fitted with a thermal Comptonization model with an electron temperature of 
50 keV and Thomson optical depth Tj ~ 1 in a slab geometry 11301 . The RXTE data 
analysis of the first 12 days of the outburst revealed a power law plus blackbody 
(with kT ~ 1 keV) interpreted as emission from the hot spot on the NS surface. Fur- 
ther simultaneous RXTE and C 'hand Va/HETGS (High Energy Transmission Grating 
Spectrometer) spectral observations collected towards the end of the 2004 outburst 
revealed a power law plus a cold black body at 0.42 keV, interpreted as the cold 
accretion disk, with no signature of the hot spot blackbody [ 134|. This is consistent 
with the lack of detected pulsations towards the end of the outburst 11291 . 

A tentative optical counterpart was identified during the 2004 outburst 01351 . 
with an R magnitude of 17.4. Follow up spectral measurements with the 4.2-m 
William Herschel Telescope on La Palma strengthened the association with the 
identification of a weak Hell and H a line 1 1361 . An optical/NIR photometric study 
in quiescence was then performed in 2005 with the 3.6-m Telescopio Nazionale 
Galileo (TNG) 033 ■ The VRIJH counterparts of IGR J0029 1+5934 were detected, 
with a strong upper limit in the K-band. The optical light curve shows variability 
consistent with the 2.46 hr orbital period. The radio to X-ray Spectral Energy Distri- 
butions (SEDs) revealed a blue component indicative of an irradiated disc 11311 . The 
SED contained also a transient NIR excess similar to that found in SAX J 1808,4- 
3658, XTE J0929-314 and XTE J1814-338. An H a line was observed in the op- 
tical spectrum of the 2004 and 2008 outbursts. Evidence was also found for an I 
band excess during quiescence in an observation performed with the 4.2-m William 
Herschel Telescope (WHT) on 2006 September 13 and 14 lfP32l . 

Radio observations taken on 2004 December 4 with the Ryle Telescope in Cam- 
bridge at 15 GHz resulted in a detection at 1.1 mJy 01381 . Follow up observa- 
tions taken one day later gave a non detection, with upper limits of 0.6 mJy at 
15 GHz Ifl39l . The Westerbork Synthesis Radio Telescope (WSRT) observed IGR 
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J0029 1+5934 between 2004 December 6 and 7, and a radio counterpart was seen 
at 5 GHz with a flux of 0.250 ± 0.035 mJy |fl39ll . Very Large Array (VLA) obser- 
vations performed on 2004 December 9 at 4.86 GHz confirmed the reported radio 
detection, with a flux of 0. 17 ± 0.05 mJy. The source was not detected in radio dur- 
ing the 2008 outburst, with upper limits of 0.16 mJy at 5 GHz in an observation 
taken between August 15 and 16 II 1401 . 



3.7 HETE J1900.1-2455 

The High Energy Transient Explorer (HETE-2) discovered HETE J1900. 1-2455 on 
2005 June 14 via an X-ray burst 01411 . Follow up observations taken on June 16 by 
RXTE/PCA showed pulsations at 377 Hz 111421 . A source distance ^4 — 5 kpc (as- 
suming a typical NS mass of 1.4M and a He rich NS atmosphere) was determined 
using the first type I X-ray burst detected with HETE-2 01431 11441. which showed 
PRE. The orbital period was determined to be 1 .4 hr with a minimum companion 
mass of O.O16M |[35]]. Unlike the other AMXPs with brown dwarf companions dis- 
cussed so far, the donor star in HETE J1900. 1-2455 is most likely a brown dwarf 
which does not require X-ray irradiation or a non-standard evolution to fill its Roche 
lobe. Simultaneous RXTE and INTEGRAL observations showed that the 2-300 keV 
spectrum is well fit by a blackbody component peaking at 0.8 keV plus thermal 
Comptonization with electron temperature of 30 keV and optical depth Zj ~ 2 for 
a slab geometry 111451 . The source was detected up to energies of 200 keV and 
the 0.1-200 keV luminosity was inferred to be 5 x 10 36 ergs~'. Although the peak 
luminosity is not particularly high compared to other AMXPs, the average mass ac- 
cretion rate (M) of HETE J1900. 1-2455 is among the highest observed due to its 
prolonged X-ray outburst that has lasted for seven years (and is still active in June 
2012) |fT46l[T47l This makes HETE J1900.1-2455 a so called "quasi-persistent 
source", i.e., an X-ray binary with a prolonged outburst duration of several years. 
Despite the long baseline provided by the outburst, pulsations were detected only in 
the first two months of the observations (Section r2.lt . 

An optical counterpart was identified with the Robotic Palomar 60-inch Tele- 
scope on June 18 J73). The source had a brightness R ~ 18.4. Near IR observa- 
tions taken with the 1.3 m robotic PAIRITEL telescope on Mt. Hopkins indicated 
a counterpart with J = 17.6 ±0.2, with no variability between observations taken 
on June 16 and 18 01481 . The counterpart was observed again with the FLWO 
1.2-m telescope at Mt. Hopkins with V = 18.09 ±0.03, R = 18.02 ±0.03 and 
/ = 17.88 ±0.02 111491 . A spectral analysis of the optical counterpart showed a Hell 
4,686 A, suggesting correct identification of the binary's optical counterpart. Fur- 
ther optical observations taken with the 8.4-m Large Binocular Telescope (LBT) on 
Mt. Graham, Arizona, showed a brightness of R = 18.64±0.02 and 18.54±0.02 on 
2007 June 13, after an initial decline of 1 .7 mag in optical, and a non-detection in X- 
rays which occurred two weeks prior I1150JI1511 . Phase resolved optical photometry 
and spectroscopy obtained in September 2006 at the 2.3-m ATT, 3.5-m WIYN and 



16 



A. Patruno, A. L. Watts 



8.4-m LBT B152I , revealed a sinusoidal modulation in the R-band, close to the or- 
bital period. However, further analysis showed that the modulation is probably due 
to a super-hump of a precessing accretion disk, and not irradiation of the donor star. 
A broad H a emission line in a 10-m Keck I spectrum taken by the same authors con- 
firmed this interpretation. No radio counterpart was identified in VLA observations 
on 2005 June 19 and 24 at 8.46 GHz lfl53l . 



3.8 Swift J 17 56.9-2508 

Swift J1756. 9-2508 was discovered by the Burst Alert Telescope (BAT) aboard 
Swift on June 7, 2007 |36l . High time resolution data from the RXTE/PCA, showed 
that the source was a 182 Hz AMXP II 1 541 . This was the fourth ultra-compact 
AMXP discovered, with an orbital period of 54.7 min. The minimum donor star 
mass, given by the X-ray mass function, is 0.0067 Mq. It has been proposed that 
Swift J1756. 9-2508 harbors a He rich white dwarf which is irradiated by X-ray flux 
from the accretor ||36l . The spectrum of the source, taken with Swift/XRT, BAT 
and RXTE/PCA, is well fit by a single absorbed power law with spectral index 
r ~ 2 |36| . The 2007 outburst length was 13 days, with no prior outbursts seen 
in archival data. A second outburst was detected on 2009 July 8 by Swift/BAT and 
RXTE/PCA 11551 . This second outburst showed a very similar length and shape as 
the 2007 outburst. No thermonuclear bursts were observed. There is evidence for a 
fluorescent Fe line in the X-ray spectrum, although the poor spectral resolution of 
RXTE prevented to constrain the physical parameters of the system 11371 . A possible 
NIR K s counterpart was identified HI 561 by inspecting the variability of 70 objects 
within the 3.5 arcsec error circle of the Swift/XRT position. The candidate counter- 
part exhibited variability of about 1.3 mag in two observations taken close to the 
end of the outburst (2007 June 19) and during quiescence (2007 July 1). No radio 
counterpart has been identified 111571 11581 and Chandra observations taken on 2007 
July 7 detected one single photon in a 1 1 ks long observations 11591 . 



3.9 Aql X-l 

Aql X- 1 is perhaps the most peculiar AMXP since it has showed pulsations at 550.27 
Hz in just 120 s of X-ray data |38|: the short baseline means that no coherent timing 
solution can be obtained. As a consequence, no orbital modulation or mass function 
has yet been measured with X-rays. The source shows thermonuclear bursts with 
burst oscillations whose asymptotic frequency is w 0.5 Hz smaller than the coher- 
ent accretion powered pulsations observed in 1998 [ 160 1 . The recurrence time of 
the outbursts is about 1 yr and they have been extensively monitored at different 
wavelengths. Twin kHz QPOs were observed at a frequency of ~ 800 and ~ 1080 
Hz 11611 . so that their frequency separation is consistent with being equal to half 
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Fig. 4 Mass-Radius relation 
(solid line) for the compan- 
ions of the four ultra-compact 
AMXPs, as determined from 
the mass function. Open 
and filled circles indicate 
an inclination of 60° and 
30°, respectively. Dotted 
and dashed curves are non- 
irradiated/irradiated white 
dwarf models, respectively. 
The numbers next to the 
curves indicate the degener- 
acy parameter. The solid thin 
curves are pure C white dwarf 
models and the numbers refer 
to the log of the temperature 
(in Kelvin). (Figure from 1 36 1) 



the NS spin frequency. The optical counterpart of Aql X-l, named VI 333 Aql, was 
identified in 1978 [ 162] as a star with B ^ 11 mag (in outburst) and B > 20 mag (in 
quiescence). The spectral type of the source was identified as a K6-M0 star II163L 
An orbital modulation at a period of 18.95 hr was reported from an optical modu- 
lation of the counterpart of Aql X-l 111641 |39l . Despite several attempts, the radio 
counterpart has been observed in only a few outbursts at a level of ~ 0. 1 — 0.5 mJy 
B1651 11661 1167L In November 2009, radio observations provided evidence for the 
emission to arise from steady jets triggered at state transitions from the soft to hard 
state Ifl68l . 



3.10 SAX J 1748.9-2021 

SAX J1748.9-2021 was discovered with the BeppoSAX/WFC on August 22, 1998 
in the globular cluster NGC 6440 111691 11701 . The distance of the source is well 
constrained thanks to the knowledge of the globular cluster distance of 8.5 ±0.4 
kpc 11711 . RXTE/PCA and BeppoSAX/WFC monitoring revealed no pulsations, but 
thermonuclear bursts were observed with no burst oscillations. The position of the 
source was coincident with the position of MX 1746-20 detected by Uhuru and 
OSO-1 in 1971-1972. However, given the large error box of these missions, it is not 
possible to establish a secure association between the two sources. SAX J1748.9- 
2021 showed a typical average X-ray spectrum with a power-law with photon index 
1.5 and a soft blackbody component at kT = 0.9 keV 1 169]. However, other spectral 
models also gave a good fit to the data 11691 . 

In August 2001, after 3 years of quiescence, SAX J1748. 9-2021 started a new 
outburst observed with the RXTEI ASM. Chandra observations allowed the estab- 
lishment of a precise position in the globular cluster and confirmed that the 1998 
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outburst was from the same source. Thermonuclear bursts without burst oscilla- 
tions were again observed. The beginning of the third outburst was detected in 2005 
March 7 during a routinely RXTE/PCA Galactic Bulge Scan Survey. The source was 
monitored until 2005 July 21 and several pointed RXTE/PCA observations were per- 
formed. The first pulsations reported for this source were seen in the 2005 outburst 
II 1721 11731 . A coherent signal at ~ 442 Hz was detected during a flux decay rem- 
iniscent of the tail of a superburst. However, neither the temporal profile nor the 
energetics of the tail were consistent with a super-burst and this lead to the sug- 
gestion that this was a new intermittent AMXP 111731 . An independent analysis of 
RXTE/PCA archival data from the 2001 and 2005 outburst revealed intermittent co- 
herent X-ray pulsations at 442 Hz, appearing and disappearing on timescales of few 
hundred seconds BT1 |40l . A further outburst was observed in December 2009 with 
RXTE and intermittent pulsations were observed again [63 1 . Observations in quies- 
cence carried out in July 2000 and June 2003 revealed an X-ray spectrum with a 
thermal component detected in both observations. The 2000 observation, however, 
required also a power law component possibly because of residual accretion on the 
NS lfT74l. 

Thanks to ROSAT/HRl archival observations and a 3.5 m New Technology Tele- 
scope (NTT) observation on 1998 August 26-27, two candidate optical counterparts, 
dubbed VI and V2, were identified II 1 701 . V2, with B ~ 22.7 mag, has now been rec- 
ognized as the true counterpart thanks to the precise astrometric position obtained 
with the 2001 Chandra observations 01751 . VLA radio observations provided strin- 
gent upper limits during the 2009 outburst and implied that the radio emission is 
quenched at high X-ray luminosities [ 176|. 



3.11 NGC6440X-2 

On 2009 June 28, the X-ray transient NGC 6440 X-2 was discovered ifTTTl in 
a Chandra observation of the globular cluster NGC 6440 (the same cluster that 
harbors the AMXP SAX J1748.9-2021). Coherent pulsations at 206 Hz were de- 
tected in a subsequent outburst in RXTE/PCA observations on 2009 August 30 [42 1 . 
The orbital period is 0.95 hr, making this the fifth ultra-compact AMXP discov- 
ered. One of the most striking features of this pulsar is the very low long-term av- 
erage mass accretion rate (M) < 2 x 10 _12 M Q yr _1 , which could be determined 
thanks to the well-defined distance to the globular cluster [69|. Peak flux is one 
of the faintest observed from the AMXPs (L2_iokeV < 1.5 x 10 36 ergs _1 ) and the 
outburst length is extremely short, with an duration of w 3-5 days. It should be 
noted that the RXTE/PCA Galactic Bulge Scan flux limit corresponds to a luminos- 
ity ~ 10 35 ergs _1 ll42l so that the presence of a longer and very faint outburst cannot 
be ruled out, although archival Chandra observations did not reveal any counterpart 
in quiescence with an upper limit on the X-ray luminosity of Lx < 10 31 ergs~'. This 
raises the question of how many faint sources of this kind might have been missed, 
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particularly important given that AMXPs seem to be associated with faint LMXBs, 
and systems of this type may therefore contain undiscovered AMXPs. 

The recurrence time of the outburst is also the shortest among AMXPs, ~ one 
month. An RXTE/PCA and Swift/XRT monitoring campaign revealed the occur- 
rence of several outbursts in 2009 (May 29-June 4, Aug. 29-Sept. 1, Oct 1-3, Oct 
28-31; I1691D . In November 2009, visibility constraints and a new outburst of the 
other AMXP SAX J1748.9-2021 prevented further observations of NGC 6440 X-2. 
In 2010 two other outbursts were detected on Mar 19 HI 781 and on June 1 1 H179II . 
In this latter outburst a strong (~ 50%rms) 1 Hz modulation was observed in the 
light-curve, which strongly resembled the 1 Hz modulation seen in SAX J 1808. 4- 
3658 EH EO). 

Despite the precise astrometric X-ray position of the source, the optical counter- 
part of NGC 6440 X-2 has not yet been identified, with B > 22 and V > 21 from 
archival Hubble Space Telescope (HST) imaging of the globular cluster when NGC 
6440 X-2 was in quiescence. During the 2009 August outburst, Gemini-South ob- 
servations did not reveal any counterpart with g' > 22. Observations carried with the 
CTIO 4-m telescope during the 2009 July outburst also did not reveal any counter- 
part with J > 18.5 andK> 17 J69). 



3.12 IGRJ17511-3057 

The X-ray binary IGR J1751 1-3057 was discovered during Galactic Bulge Moni- 
toring by the INTEGRAL satellite in 2009 September 12 111801 . Coherent pulsations 
at 245 Hz and an orbital modulation of 3.5 hr were measured with follow up ob- 
servations by the RXTE/PCA [43|. Given the X-ray mass function, the minimum 
donor star mass is O.13M . The X-ray light-curve was also extensively monitored 
by Swift/XRT, XMM-Newton, RXTE/PCA and Chandra. The light-curve showed a 
typical fast rise with a slow decay lasting for about 20 days, before entering a fast de- 
cay phase [181, 1 82 1 . Curiously, while the source was in the fast decay phase, a very 
short and faint outburst of the ultra-compact AMXP XTE J 175 1-305 appeared in 
the same field of view (FoV) of RXTE. The mini-outburst of XTE J 175 1-305 lasted 
less than 3 days, and pulsations were clearly visible at 435 Hz 111 101 . Thermonuclear 
bursts were observed in IGR J1751 1-3057 with Swift 111831 . and burst oscillations 
were detected in RXTE data 1 184, 11821 . None of the bursts showed PRE, so there 
are several different distance determinations in the literature. The most constraining 
gives an upper limit on distance of < 5 kpc II 1 821 . 

Swift spectral analysis revealed a power law (F ~ 2) plus a blackbody at 0.9 
keV lfT8l . Broadband spectral data from simultaneous RXTE/PCA and Swift/XRT 
observations could be well fitted with a three model component 11 1 851 : a disk with 
low temperature (kT ~ 0.24 ke V), a hot black body (hot spot with kT ~ 1 ke V) and a 
Comptonization component originating from the accretion shock (electron temper- 
ature T e ~ 30 keV; Zj ~ 2). The spectral fitting also required a fluorescent iron line 
at 6.4 keV with Compton reflection and provided an interstellar absorption column 
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of Nh = x 10 22 cm- 2 . Very similar results were obtained using simulta- 

neous XMM-Newton and RXTE/PCA spectral data |fT8ll . Combined RXTE/PCA, 
Swift/XRT and INTEGRAL data also gave consistent results using a thermal Comp- 
tonization model II 1 861 . Twin kHz QPOs were observed at the beginning and at the 
end of the outburst, and had Av ~ 120 Hz, about half the value of the spin fre- 
quency II 1 8711 . A NIR counterpart with brightness K s — 18.0 ±0.1 was identified 
with the 6.5 m Magellan Baade telescope on 2009 September 22 [188]. A second 
observation on October 7, while the X-ray flux was fading rapidly, showed no coun- 
terpart with 3 a upper limits of K s > 18.8. No radio counterpart has been found, with 
upper limits of 10.0 mJy [ 1 89 1 . 



3.13 Swift J 1749.4-2807 

The first observation of Swift J1749. 4-2807 was made on June 2, 2006 with the 
Swift/BAT telescope, which recorded a high level activity consistent with an uniden- 
tified source 01901 . This was later identified as a "burst-only" accreting NS binary 
(i.e., a NS that is bright enough to be detectable only during the occurrence of ther- 
monuclear explosions) thanks to a spectral analysis of the SwiftfBAT data 11911 . 
The analysis was consistent with the observation of a thermonuclear burst for a 
source distance of 6.7 ± 1.3 kpc. Further Swift/XRT data analysis revealed an X- 
ray counterpart of the burst and archival XMM-Newton data showed a faint point 
source coincident with the Swift/XRT position 11921 |19ll . Increased activity was 
reported in 2010 by INTEGRAL, and was linked to the accretion powered emission 
process 11 1 9 3 1 1 1 94l . Pulsations were discovered immediately afterwards at 518 Hz in 
RXTE/PCA follow-up observations I1195II451 . A very strong second harmonic was 
detected at 1036 Hz 1 196 1 and orbital modulation was measured at 8.82 hr, giving 
a minimum donor mass of about 0.6M^ lfl97l[T98l . Swift J1749.4-2807 was finally 
identified as the first eclipsing AMXP [199| and the inclination of the system was 
constrained to be i ~ 74° — 78° 1 199, 45]. The binary showed three eclipses with a 
duration of 2172 ± 13 s, allowing the first detection of Shapiro delay effects in X-ray 
timing data for an object outside the Solar System 11 1 991 . 

To fit the 0.5-40 keV spectrum, an absorption column density of Nh = 3.0 x 
10 22 cm~ 2 and a power-law with spectral index F ~ 1.7 are required [200] . The 
absorption column is about 3 times larger than the expected Galactic column density 
in the direction of the source. The X-ray light-curve showed an exponential decay 
with a Swift/XRT non detection after 11 days since its first 2010 observation. Due 
to the crowded field (the source is on the Galactic plane), no single counterpart has 
been identified in NIR counterpart searches [70|. More than forty counterparts were 
identified in the X-ray error circle with ESO's Very Large Telescope (VLT). 
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The most recent addition to the AMXP family was discovered by INTEGRAL on 
2011 August 11 B201I . Soon after, coherent X-ray pulsations were discovered in 
RXTE/PCA data at 401 Hz with an orbital period modulation of 3.8 hr [46]. Type 
I X-ray bursts were observed with INTEGRAL [202 1 and RXTE. Burst oscillations 
were detected [203 204] and evidence of PRE placed the source at a distance of 
7.6 kpc [ 203 1 . Swift/XKT observed the source returning to quiescence on 2011 
September 19 12051 . A candidate NIR counterpart of IGR J 17498-2921 was first 
detected in archival data taken with the 4-m VISTA telescope at the Paranal ob- 
servatory [206]. The optical counterpart was detected on 2011 August 25 and 26 
with the 2-m Faulkes Telescope South [207 1 . The position detected was consistent 
with both the Chandra and the NIR counterparts. However, further observations 
suggested that the NIR/optical counterpart is a foreground star aligned by chance 
with the X-ray source [208]. This counterpart was confirmed to be a foreground 
source with observations taken at the 2.5-m Irenee du Pont telescope. The analysis 
suggested that, given the distance and extinction of IGR J 17498-2921, the expected 
magnitudes are R ~ 28 and J ~ 23.3, both well above the limit reached by the opti- 
cal/NIR observations 12091 . 



4 Accretion Torques 

Once the donor star overflows its Roche lobe, gas flowing through the inner La- 
grangian point L\ carries large specific angular momentum and hence forms an ac- 
cretion disk around the NS. The type of disk depends on the microphysical condi- 
tions governing the gas dynamics 12101 121 ll 12121 12T31 . If gas pressure dominates, 
the disk will be geometrically thin and optically thick 12141 121511210 1 with material 
moving in Keplerian orbits with orbital frequency: 

1 [GM wtj / M \ 1/2 f R V 3/2 

At a distance of a few tens of kilometers from the NS, the gas orbits several hundred 
cycles per second and flows almost undisturbed until the magnetic field of the NS 
(that for AMXPs is of the order of 10 8 G) is strong enough to perturb its orbit. At 
the magnetospheric-radius r m , the kinetic energy of the free-falling gas becomes 
comparable to the magnetic energy of the NS magnetosphere : 

where M is the mass accretion rate at the inner disk boundary, /I is the dipole mag- 
netic moment of the NS and r& is the Alfven radius. This latter parameter is calcu- 
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lated assuming spherical accretion: 

^\ piM . (3) 

The parameters p and v are gas density and velocity, respectively. The term % w 
0.3 — 1.0 is a correction factor due to the non-spherical geometry of the problem 
and is required because the gas orbits in a disk rather than falling radially from 
every direction. In the disk geometry, magnetic and fluid stresses balance when: 

B p B$r 2 =M^y- (4) 
or 

where B p and B^ are the poloidal and toroidal components of the NS magnetic 
filed, r the radial coordinate measured from the NS center and the azimuthal 
velocity of the plasma at r. If one assumes that the transition region Ar connecting 
the unperturbed plasma flow far from the NS and the magnetospheric flow is much 
smaller than r m , then the above equation takes the form: 

B p B^r 2 m Ar^Mv^r m (5) 

Once the gas reaches the transition region Ar, it stops flowing in Keplerian orbits 
and starts to co-rotate with the magnetosphere. The gas exchanges angular momen- 
tum with the magnetosphere and changes the NS spin frequency. The NS is spun 
up if its specific angular momentum is smaller than that of the accreting gas, and 
otherwise spun down. The spin-up/spin-down condition can be thought of in terms 
of characteristic radii: if r m is smaller than the radius where the Keplerian frequency 
equals the NS spin frequency (the co-rotation radius r co ), the NS is spun up'; other- 
wise it is spun down. The co-rotation radius can be defined as: 

r co = 1 683Mj/ 3 v7 2/3 km. (6) 

It is important to stress that this is an over-simplified picture of the true physical 
conditions close in the inner disk. In this simple description, the accretion torque 
exerted on the NS for a Keplerian disk truncated at r m , with r m < r co , is: 



N = 2kIv s = M\jGMr m (7) 

where / is the moment of inertia of the NS, V s the NS spin frequency derivative and 
G the universal gravitational constant. At radii larger than the co-rotation radius, 
the magnetic field lines are threaded into the accretion disk and dragged by the 
high conductivity plasma so that an extra torque due to magnetic stresses has to be 
expected 121 1II2121|216||2T71 in addition to the torques due to the matter flow. Spin 
down due to dipole emission is also always present. A possible way to express the 
total torque acting on the pulsar is 02181 : 
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where Edipole lS the energy loss due to dipole radiation and n((0) ~ ± 1 is a dimen- 
sionless function that depends on the fastness parameter co = (r m /r t -„) 3 ^ 2 : 

n(co) = taiihfi^J (9) 

This term takes account of the gradual transition from the spin-down to spin-up zone 
in the accretion disk 12171 12181 . The two extra terms in Eq.® have a minor effect 
during most of the outburst, when the mass accretion has the largest weight in deter- 
mining the net torque. The expression in Eq.(|7) is therefore a good approximation 
most of the time. 

As the AMXP is spun up, r co moves towards and eventually reaches r,„. When 
this happens, the pulsar is said to have reached the "equilibrium spin period" P ecl : 

— 3/7 

P eq * BfM-J 1 ( —JL—^ R\W ms (10) 
Substituting into Eq.dTOl the surface magnetic field derived from dipole spin down 112191 



B — \l — * 3.2 x IO^gVpJm^^- (11) 

w 8k 2 R 6 sina 14 sina 



(where we have assumed / = 10 45 gcm 2 , R = 10 km and a is the misalignment 
angle between spin and magnetic axes), one obtains a relation between P eq and P. 
When the accretion rate reaches the maximum Eddington rate, this relation defines a 
"spin-up line" in the P — P diagram of radio pulsars, above which millisecond radio 
pulsars should not be found. Indeed, if millisecond pulsars are created in LMXBs via 
accretion torques, then the maximum possible torque is set by the Eddington limit. 
As of today, the spin-up line has never been crossed by any radio pulsar, with the 
exception of the y and radio pulsar PSR 1823-3021A 112201 in the globular cluster 
NGC 662^1 (see FigureEJ. 

To detect the effect of accretion torques in AMXPs one simply needs a mea- 
surement of V s , which in principle is straightforward when using coherent timing 
techniques. Unfortunately, however, determining the location of the magnetospheric 
radius is a non-trivial problem. 

The main reason is that Ar/r m is assumed to be much smaller than unity, 
which is a good approximation only for slow accreting pulsars with high 



2 Note, however, that the spin-up line depends on several parameters which are difficult to constrain 
like the angle a. Its position in the P — P diagram is therefore subject to uncertainties (see 112211 
for a discussion). 
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Fig. 5 P — P diagram of radio pulsars with a spin up line (dashed black line) above which radio 
pulsars should not be found if they are "recycled" via accretion. The pulsar PSR 1823-3021A 1 220 1 
is not shown in the figure. The red lines identify magnetic fields of different strengths. Data from 
the PSRCAT on-line data archive. 



magnetic fields. This might not be the case for AMXPs with rapid rotation 
and weak magnetic fields (see for example Eq. 5 and the underlying assump- 
tions). 



Furthermore we have considered so far only gas-pressure dominated accretion 
disks, considered to be a good assumption at the low average accretion rates of 
AMXPs. However, when the outbursts are close to their peak luminosities, radiation 
pressure may play a role in the exchange of angular momentum between the gas 
and the NS. At these accretion rates accreting plasma and magnetic field lines may 
couple and modify the amount of enhanced angular momentum B222I [223], thus 
invalidating the use of Eq.((7][8]l. 

Noting these caveats, theoretical expectations for the average spin frequency 
derivative of an AMXP can be calculated with the simplest accretion model (Eq. |7]i 
once one has a reasonable estimate of the mean mass accretion rate M during an 
outburst. As discussed, such a measure is difficult to obtain, but several estimates of 
M are present in the literature which can be taken as a first step to compare the mea- 
sured v v with the expected values (see for example [224 1 and [223]). The explicit 
expression for the expected vf v can be obtained by substituting Eq.© into Eq.©. 
The substitution gives: 
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vfP = 2.3 x 10- x % x l 2 M%Ml n A BfR^Uzs- x (12) 

As one can easily verify, the calculated values of vf* are all of the order 
of lO-^Hzs -1 when neglecting the fact that at low mass accretion rates, when 
r m > f co the "propeller regime" can set in, with the spin-down terms dominating 
in Eq.®. In this phase, the specific angular momentum of the accreting plasma is 
insufficient for spin up to occur, and the centrifugal barrier removes angular mo- 
mentum from the AMXP, slowing it down. It was originally suggested centrifugal 
inhibition by the rotating magnetosphere would expel gas from the system and shut 
down the accretion process 112251 . In fact for this to happen, r m must exceed r co by 
a margin of at least 1.3 so that matter can be accelerated to the escape velocity and 
can be flung out of the disk |[84ll217|[85l . Accretion can in fact still take place, and 
two propeller regimes have been observed in 3D-MHD simulations of accreting pul- 
sars 1 226 227 1 : a "weak propeller" with no outflows and a "strong propeller" with 
expulsion of material. In both cases channeled accretion is still ongoing so that a 
spinning down AMXP could be in principle observed via coherent timing measure- 
ments in either propeller regime. The propeller is therefore expected to affect V s , 
but determining its onset is a very difficult task. Indeed, AMXP observations do not 
help much to constrain the theory in this case, since the only evidence for possible 
propeller phases comes from SAX J1808.4-3658 and IGR J00291+5934. In the first 
source a propeller may explain a strong 1 Hz modulation observed towards the end 
of the outburst ([80], Section lXTT i and pulse phase lags ( 12281 . Section lBTTT i whereas 
for IGR J0029 1+5934 the propeller model may explain the anomalous double 2008 
outburst ([33], Section l3~6l l. In both cases further evidence is required to confirm 
a propeller phase, and there is no other conclusive observational evidence that any 
AMXP enters this phase. Coherent timing does not provide further insights, since 
the propeller is expected to start when the mass accretion (and thus the source lumi- 
nosity) drops substantially, with a consequent decrease of the S/N of the pulsations. 
Furthermore the propeller may not last sufficiently long to allow measurement of 
the spin frequency derivative. An interesting claim that the onset of the propeller 
had been observed in the 2002 coherent timing data of SAX J1808. 4-3658 was pro- 
posed [229 1 to explain a pulse phase jump that occurred at the onset of the fast decay 
of the outburst. However, the pulse phase jump has been interpreted differently in 
other works lfl4l 1911 1401 as it might be severely contaminated by X-ray timing noise. 



4.1 Coherent Timing Technique 

Coherent timing analysis of the phase evolution of AMXPs can be performed after 
converting photon times of arrival from the spacecraft non-inertial reference frame 
to the Solar System barycenter (an approximate inertial reference frame) and cor- 
recting for general and special relativistic perturbations due to planets and minor 
bodies in the Solar System that affect photon propagation. Before discussing obser- 
vational tests of accretion theory it is useful to clarify the observables that play a 
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role in coherent timing studies. AMXPs show pulsations that are too weak to detect 
as single pulses: to reconstruct the signal it is necessary to fold the data in segments 
of several hundred seconds to obtain a pulse profile that is the average of several 
hundred thousand NS cycles. It is then possible to measure the fractional ampli- 
tude of the pulsations, and the pulse phase. Pulse profiles of AMXPs are generally 
highly sinusoidal, with little or no harmonic content beyond the fundamental fre- 
quency. In this case the pulse phases can be measured by choosing a fiducial point 
(e.g. the pulse peak) and tracking the variation of the phase at this point over time. 
Sometimes, however, strong harmonic content is observed, with the pulse profile 
shape varying during an outburst. This means that unlike in radio pulsar pulsar tim- 
ing, where average profiles are often very stable, there is no stable fiducial point 
on which to base timing analysis. To avoid this problem it has become standard to 
decompose the pulse profile into its harmonic components and measure pulse am- 
plitude bfr and phase fa of the k-th harmonic (k = 1 for the fundamental, k = 2 for 
the second harmonic and so on) via the expression: 



where bo is the unpulsed component, xj is the number of counts detected in the j- 
th bin of the pulse profile, with j — 1,2, ...,N. This way the fiducial point of each 
harmonic is well defined since all harmonic components are pure sinusoids. The 
fractional amplitude can be measured by adding in quadrature the fractional ampli- 
tudes of each harmonic: 



where is the fractional amplitude of the k-th harmonic, N p t, — Y.j x j the total 
number of photons in a profile and B the total number of background photons. Note 
that this definition gives a fractional amplitude larger by a factor y/2 than the of- 
ten reported rms fractional amplitude. We use this definition of fractional amplitude 
since it has an immediate physical meaning as the pulsed flux fraction. To avoid 
confusion, we always refer to the fractional amplitude as "sinusoidal fractional am- 
plitude" when we use the definition given in Eq.dT4b or otherwise to "rms fractional 
amplitude". Note that in the AMXP literature "rms fractional amplitude", "sinu- 
soidal fractional amplitude" and "peak-to-peak fractional amplitude" are also used. 
This latter is calculated by measuring the flux at the peak and at the minimum of the 
pulse and dividing by the average flux. 

Pulse phase can be measured as a function of time and each harmonic analyzed 
separately. We drop the subscript k from the phase symbols, since the following 
equations are valid for each harmonic. The pulse phase series can be decomposed 
into several term encoding different physical effects: 




(13) 




(14) 



0) = 0L + 02 + Qorb + 04 + §M + 0iV 



(15) 
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<j>L = <j>o + vt is the linear term due to a constant spin of the NS (0o is an initial 
reference phase) 

= \ vt 2 is the quadratic variation due to constant spin up or spin down 
(j) or j, contains the effect of the orbital motion of the NS around the companion 
and, as a first order approximation, can be calculated by measuring the delays in 
the time of arrival of photons: 



Ai . . . 

tern = tarr Sim Sin 

C 



2% 

\Jarr~ *asc) 



(16) 



where A i is the semi-major axis of the NS orbit, c the speed of light, i the inclina- 
tion of the orbit with respect to the observer, T asc the time of passage through the 
ascending node, Pi, the orbital period and t em and t urr the photon emission and ar- 
rival time 1 230 1 . In this expression we have assumed that the orbit is perfectly cir- 
cular, a good first order approximation in all AMXPs. More sophisticated models 
for nearly circular orbits can be found for example in [123 11 . 

• §a gives phase variations related to uncertainty in astrometric source position, 
which introduces a spurious frequency and frequency derivative offset. These 
offsets can be expressed as dlMl t: 

Av = Voe(a e cosj5/c)(27z/P e )cosT (17) 
Av = -Voe(fl©cosj3/c)(2^// , ®) 2 sinT (18) 

Here Vo is the true pulse frequency, £ the position error parallel to the plane of 
the ecliptic, j3 the ecliptic latitude of the AMXP, a^ and P® the Earth semi-major 
axis and orbital period and T = Int/P^ the orbital phase of the Earth. Phase zero 
is defined as the point where the Earth is closest to the AMXP and for order of 
magnitude estimates can be assumed such that cost = sinr = 1. 

• 0m refers to unavoidable phase wandering due to measurement errors, and is 
normally distributed with an amplitude predictable by propagating the Poisson 
uncertainties due to counting statistics. 

• 07v is a subtle term covering residual phase variations that do not fall into any 
of the previous categories. These are usually called "X-ray timing noise", by 
analogy with the timing noise often observed in radio pulsars. Note that if the 
spin-up (or spin-down) process is not constant in time we do expect terms higher 
than the quadratic ((jig) and these are considered, in our definition, as part of 0# 
even if they are true variations of the NS rotation. This should be expected if, for 
example, the accretion torques exhibit stochastic variations. 

When performing coherent timing after having removed the <p or t, term and as- 
suming that fa is negligible, the observer measures a pulse frequency and its time 
derivatives. These can be determined using a Taylor expansion: 

Ht) = <h + ^ it - to)+ d ^t^t + ... (19) 
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The parameters V = V = etc., can be determined by fitting the pulse phases 
with standard x 2 minimization techniques. 

The observable quantities here are not the spin frequency v s and derivatives, 
but the pulse frequency v and derivatives which are encoded as a combina- 
tion of fa, fa and fa. The pulse frequency is the frequency of the pulsations 
detected by the distant observer, whereas the spin frequency is the rotational 
rate of the NS as measured by the distant observer. The assumption that pulse 
and spin frequency (and derivatives) are identical may not always be true. For 
pulse and spin frequencies to be identical, fa must have no linear compo- 
nent. Similarly for the pulse and spin frequency derivative: only if fa has no 
quadratic component will the two be the same. 



4.2 Observations: Accretion Torques in AMXPs 

With the exception of the intermittent pulsar Aql X-l, all of the AMXPs have shown 
pulsations of sufficient quality and with a sufficiently long baseline to constrain the 
pulse frequency derivatives. As a rule of thumb, the condition that must be met 
to detect a pulse frequency derivative v in a data segment of length At is that 
Cms < v(4f) 2 , where (J rms is the root-mean-square error on pulse phases. Since 
all AMXPs (bar Aql X-l) have a rms ~ 0.01 cycles, and baselines of several days, 
we are sensitive to v ~ 10~ 15 — 10~ 13 Hzs _1 . These values overlap the range of 
expected V given by Eq.© for typical M and weak dipolar B fields ~ 10 8 — 10 9 G. 

Several papers have reported pulse frequency derivatives in AMXPs. The first 
measure was made for the source IGR J0029 1+5934 [30 1 with a reported V = 8.4 x 
10 _13 Hzs _1 during its 2004 outburst. Table [3] summarizes the measurements for 
all AMXPs. However none of these values, including the measurement made for 
IGR J0029 1+5934, take into account the presence of X-ray timing noise in the pulse 
phases. One must therefore bear in mind that what is reported is the combined effect 
of fa and fa, as explained in Section 14.11 Nevertheless timing noise has different 
strength in different sources, so that pulse and spin frequency derivatives may in 
some cases not be so dissimilar. We have marked each AMXP in Table[3]with a "w" 
for weak timing noise and "s" for strong. This distinction is somehow arbitrary, but 
is useful to understand that probability of the pulse frequency derivatives differing 
from the true spin frequency derivatives. 

To verify whether the pulse frequency derivative is a robust indicator of the spin 
frequency derivative v v , one has to investigate the strength of fa on the timescales 
over which v is measured. To do this, one can use Monte Carlo (MC) simulations 
to estimate the true uncertainty on the pulse frequency derivative [14|. The MC 
method works as follows: one fits a second-order polynomial to measure v and v 
with standard % 2 minimization techniques after having removed all other effects 
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Table 3 Pulse Frequency and Spin Frequency Derivatives 



Source 


Pulse free], deriv. V 


Timing noise PFD 


Spin freq. deriv v v SFD 




(PFD) [Hzs- 1 ] 


strength Reference 


(SFD)[Hzs _1 ] reference 


SAX J1808. 4-3658 


4.4 x 10 IJ ; -7.6 x 10" 


14 s 


2291 


< |2.5| x 10- 14 


[14J115J 


XTE J 1807-294 


2.5 x 10- 14 


s 1251 


< 141 x 10- 14 


2321 


IGR J0029 1+5934 


[5; 11] x 10-" 


w 130 31 34 +5.5 x 10 13 


30]|3T][32]|33][34] 


XTEJ1814-338 


-6.7 x K)-' 4 


s 


28] 


<|1.5|xl0- 14 [233] 


XTE J0929-314 


-9.2(4) x 10 14 


s 


[21] 


XTE J 175 1-305 


5.6 x 10-' 3 


w 


m 


IGR J1751 1-3057 


1.6x 10-° 


s [44] 






IGR J 17498-2921 


-6.3 x K)-' 4 


w no 






HETEJ1900.1-2455 




s 




[2.3;0.4] x 10-" [62] 


Swift J1756.9-2508 




w 




<|3|xl0- 13 [37] 


Swift J1749.4-2807 




s 




< 1.2|xl0- 12 fl99l 



(such as orbital variations) and obtains an estimate of v ± 0" v and v ± Cy. If in the 
phase residuals there are still substantial variations of the pulse phase, in excess 
of that expected from measurement errors alone (i.e., 0m) - which can be verified 
by checking whether % 2 gives a statistically unacceptable fit - then the statistical 
errors 0" v and o~y are not good representations of the true uncertainties of the spin 
parameters. Instead one can take the phase residuals, calculate a power spectrum 
and simulate several thousand time series with nearly identical noise content as the 
original phase residual time series lfl4l . At this point v and v can be measured 
for each simulated time series, and a distribution of parameters constructed. The 
standard deviation of the distribution of v and v provides a good representation of 
the true uncertainties on the spin parameters. In this way one can immediately check 
whether the noise content fa affects the measured value of v and v. Applying this 
technique has revealed discrepancies between pulse and spin parameters in several 
AMXPs (Table [3]). The most striking finding has been the non-detection of spin 
frequency derivatives in several AMXPs, with upper limits below the expected Vj. 

The use of MC simulations is not the final word on this problem, since the ori- 
gin of timing noise remains unexplained. It also has its limitations: if the lowest 
Fourier frequency component of timing noise is comparable to the length of the 
data segment over which v and v are measured, MC simulations cannot distinguish 
pulse and spin parameters. It has been noted that the phase residuals obtained af- 
ter removing v and v from the pulse phases of two AMXPs (XTE J18 14-338 and 
XTE 1807-294) are anti-correlated with variations in X-ray flux Il28ll25l . The anti- 
correlation improves substantially if one fits a simple v = const model, suggesting 
that timing noise, which is related to the X-ray flux variations, is almost entirely 
responsible for the measured v 123411 . Correlations (or anti-correlations) between 
pulse phase residuals (with respect to a v = const model) and X-ray flux have now 
been found in at least six AMXPs where such studies have been carried out 02341 . In 
some cases the correlation was striking (e.g., in XTE J1814-338 12331 ) leaving little 
doubt that the pulse frequency derivative v is not the spin frequency derivative v s . 
This discovery also suggests that it is the pulse phase <j) and not its second deriva- 
tive (v) that correlates with the flux (and thus the mass accretion rate), as would be 
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expected from Eq.0. One way to test these findings is to measure "instantaneous" 
short timescale spin frequencies and check whether these scale with the bolometric 
flux as Vj x Ffo p where y is a scale factor that depends on accretion disk structure. 
The bolometric flux is, however, almost never available, since observations usually 
cover only a narrow (high) energy band. If we assume that the X-ray flux is Fx x Fj, i 
then we expect to see v 4 °< F^. However this assumes that the mass accretion rate 
M oc Fx, which has been shown to be untrue in some LMXBs [235 1. There is one 
additional complication: AMXPs have outbursts that are almost always too short to 
split data segments long enough to yield more than one short-term "instantaneous" 
spin frequency derivative. 

If one keeps these caveats in mind, then it is possible to test the relation between 
V s and F x for some AMXPs, such as XTE J 1807-294 and XTE J18 14-338 which 
have shown long outbursts (~ 100 and ~ 40 days, respectively) and high S/N ratio 
for the pulsations. In XTE J 1807-294, the "instantaneous" V changes sign several 
times during the outburst decay and all "instantaneous" spin frequency derivatives v. s 
are insignificant and consistent with being part of the underlying timing noise pro- 
cess 112321 . In XTE J814-338, the instantaneous v are too large (up to 10~ n Hzs" 1 ) 
to be physically meaningful, requiring accretion rates well above the Eddington 
limit [236]. This is at odds with the fact that AMXPs are faint and rarely reach 
accretion rates above 10% Eddington, suggesting that V is different to v s and is 
completely dominated by timing noise. This does not prove that v s does not scale 
with flux, but only that measurements of "instantaneous" v s are contaminated by 
timing noise, which prevents testing of the relation between V s and Fx- 

Theoretical expectations for V s mostly exceed observed values measured by tak- 
ing into account the contamination of timing noise (Table [3). This means that most 
AMXPs do not behave in accordance with the predictions of accretion torque theory 
(Eq|8j. In addition, once gas attaches to the magnetic field lines of the AMXP, an 
exchange of angular momentum is inevitable unless r m = r co . But this cannot be 
the case throughout an outburst: r m M~ 2 / 7 and the mass accretion rate varies, so 
r m must at some point differ from r co . Larger torques, which are expected close to 
the outburst peak, are clearly ruled out in some AMXPs lfl"4l l26l 1233 1 . The angular 
momentum of the accreting plasma must however be transferred somewhere so that 
the mismatch between expected and observed v s appears problematic. 

4.2.1 The Origin of X-ray Timing Noise 

The phase wandering of pulsations in coherent timing analysis requires careful con- 
sideration if we are to understand how a NS responds to accretion. Theoretically 
the phenomenon is easy to explain if one assumes that the accretion hot spot is 
not completely anchored to one location on the surface. One possibility is that the 
magnetic and rotational axes of the NS are almost aligned and the hot spot wan- 
ders around the magnetic axis by small displacements, generating large variations 
in phase and amplitude 12381 . In this case pulsed fractions and times of arrivals 
should be anti-correlated. A moving hot spot has been observed in 3D-MHD sim- 
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ulations of magnetized accreting NSs 1 2391 12371 . In particular, it was shown that 
the usual assumption of a fixed hot spot is valid only for large misalignment angles 
between the magnetic and spin axes (Figure|6] [237 1), in good agreement with the 
magnetic and rotational axis alignment model [240 238 1. This strengthens the idea 
that most X-ray timing noise may be related to a moving hot spot configuration. 
Alternatively, the phase wandering might be related to complexity in the structure 
of the pulsar's magnetic field, with different multipole components dominating the 
accretion process at different accretion rates [12411 . 

Timing noise may also be due to fluctuating accretion torques [242, 243 1. This 
has been studied for strongly magnetized accreting pulsars like Vela X-l [ 244 1 , 
and it is natural to expect something similar in AMXPs. However as discussed, 
this model requires unphysical large accretion rate fluctuations to explain all timing 
noise or the large jumps in phase observed for example in XTE J 18 14-338 (and 
probably also SAX J1808.4-3658 and XTE J1807-294, Eg] |25] EH [2361 COD)- 
Nonetheless some X-ray timing noise, or even its entirety for AMXPs with weak 
timing noise content, could in principle be related to accretion torque fluctuations. 
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Fig. 6 Evolution of the position of the hot spot during accretion (v s = 4.1ms). [l represents the 
magnetic dipole axis. Note how the densest portion of the circular hot spot (encoded by the p 
scale, bottom right) shifts during time even if the observer is in the rotating reference frame of the 
NS (Figure from (237J.) 
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Fig. 7 Mass-Radius relation for several different EoS of ultra-dense matter, and constraints on 
the M-R relation of the AMXP SAX J 1808.4-3658. These result from modeling pulse profiles 
observed in the 1998 and 2002 outbursts (Figure from 112491 ). 

5 Pulse Profiles 

AMXP pulse profiles form in regions of strong gravity and encode information 
about the physical properties of the emitting regions and the compactness of the 
NS H2451 12461 12471 12481 12491 . Light bending, aberration and relativistic Doppler 
shifts all affect pulse profile shape [250, 251 252 , 253 1 and can be measured using 
coherent timing analysis. This can be used to constrain the EoS of ultra-dense matter 
(Figure|7]i although large uncertainties still exist: in part due to model-dependencies, 
but more importantly because of the low S/N of the relativistic features. 

Pulse profile shapes can also constrain the geometry of the NS magnetic field 
[254, 255 256 1 since multipoles generate different accretion columns and hot spots 
on the NS surface (see, however, [257], for a recent discussion of the problem). The 
information contained in pulse profiles can be extracted by studying pulse ampli- 
tudes, the harmonic content of the signal and the time of arrival of the pulsations. 
Double peaked pulse profiles, for example, may be observed when two antipodal 
hot-spots exist on the NS surface 1125 11 (238 . 253 1; a high harmonic content is sugges- 
tive of a complex field geometry [258 1 and the energy dependence of the pulse am- 
plitudes and the time lags of the pulsations provide a way to explore the angular pat- 
tern of the radiation emitted from the NS surface [245. 1 12. 253 1. The temporal evo- 
lution of the pulse shapes also reveals details of the complex disk/magnetospheric 
interaction responsible for magnetic channeling of the accreting plasma [259] [91]. 
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5.1 Pulse Fractional Amplitudes and Phase Lags 

The pulse profiles of AMXPs are sinusoidal, with only the fundamental frequency 
detected in most cases (including IGR J00291+5934, XTE J0929-3 14, XTE J175 1- 
305, SAX J1748.9-2021 and IGR J17498-2921). In some AMXPs there is a first 
overtone (e.g. XTE J18 14-338), and higher overtones are seen only sporadically in 
just three AMXPs (SAX J1808.4-3648, Swift J1749.4-2807 and XTE J1807-294). 
Pulse fractional amplitudes reach values of a few percent, with the highest sinu- 
soidal amplitude measured being 30-40% (Swift J1749.4-2807 and XTE J1807- 
294) and the lowest being 0.43% (HETE J1900.1-2455). With the exception of 
Swift J1749. 4-2807, overtones do not generally contribute to the amplitude of the 
pulsations by more than about ~ 5%. 

All AMXPs show an energy dependence of the pulse fractional amplitudes. 
In some sources the fractional amplitude increases steeply with energy, e.g. Aql 
X-l [33 Swift J1756.9-2508 |26| and SAX J1748.9-2021 00). In this latter 
source, the fractional amplitude increases by a factor ~ 10 between 2 and ~ 
20 keV. Some AMXPs have fractional amplitudes that drop with energy (SAX 
J1808.4-3658 [260|, XTE J175 1-305 \U^,XTE J0929-3 14 EH , HETE J 1 900. 1- 
2455 ED, IGR J1751 1-3057 (9TJ[H6|). Other AMXPs (IGR J00291+5934 and 
XTE J 1807-294) are more complex, with the amplitude rising and falling in differ- 
ent energy bands [30 26 1. The energy dependence of the pulsed fractions provides 
insights into the pulse formation process, since a simple hot spot emitting blackbody 
radiation with a temperature contrast with respect to the NS surface would produce 
pulsations which increase in amplitude at higher energies in the observer reference 
frame 1125 1 1 - 

Another interesting finding has been the correlation or anti -correlation between 
the fractional amplitude of the fundamental (where most of the pulse power is ob- 
served) and the time of arrival of the pulsations of some AMXPs. The most promi- 
nent example is XTE J 1807-294 where high/low fractional amplitude pulses arrive 
systematically earlier/later than predicted by the timing model. This finding, and 
a similar anti-correlation between the time of arrivals and the X-ray flux, may be 
evidence for a hot spot moving on the NS surface ||26) . Similar conclusions were 
reached by other authors to explain the timing noise observed in the timing residu- 
als of several AMXPs |[T4l |25 28 1. Some of these observations can be understood 
in terms of the hot spot wandering model described in Section 14.2.11 although no 
final confirmation of the models has yet been put forward. Despite many open prob- 
lems that need to be solved (such as why such correlations are not observed in 
AMXPs like SAX J1808.4-3658) the moving hot spot model help us to understand 
the AMXP pulse formation process. 

Soon after the first AMXP discovery, it was realized that the lowest energy 
("soft") photons that compose the pulsations arrive on average later than the high en- 
ergy ("hard") photons 112601 . In SAX J1808.4-3658, hard photon arrival times tend 
to saturate at some energy E > 10 keV, with soft photons at about 2 keV accumulat- 
ing a lag of about 0.08 rotational cycles (or 200 /is). Similar behaviour was soon dis- 
covered in other five AMXPs: XTE J0929-314 (0.14 cycles ED), XTE J175 1-305 
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(0.06 cycles, Ol), XTE J1814-338 (0.016 cycles, IB5TT ). XTE J1807-294 (0.1 
cycles, ED) and in the intermittent pulsar HETE J1900.1-2455 (0.07 cycles, |6ll ). 
In IGR J00291+5934, soft lags were observed with hard photons leading by 0.06 
cycles at 6-8 keV |3Q| . Above this energy, harder photons (E > 8 keV) start to grad- 
ually reduce their lead until at energies of 30-100 keV there are no phase lags with 
respect to the soft photons (at ~ 2 keV). The origin of these lags is poorly under- 
stood, and has been discussed in terms of the different angular distribution ("fan" 
and "pencil" beams) of a two spectral component model 112451 or in terms of Comp- 
ton down-scattering of hard X-ray photons from the cold disk plasma or the NS 
surface lfl45l . Detailed study of the soft lags in SAX J1808.4-3658 12281 suggests 
that they have a flux dependence, with the lag being almost zero at high fluxes, in- 
creasing steadily at lower fluxes and then decreasing again below a flux threshold 
coincident with the onset of the rapid decay. This may be linked with the changing 
properties of the accretion disk as it transitions towards the propeller phase. 



5.2 Pulse Shape Evolution 

Most AMXPs with strong harmonic content show pulse shape variability that often 
correlates with the stage of the outburst. The pulse profiles observed in different out- 
bursts of SAX J1808. 4-3658 have characteristic shapes that can be associated with 
the specific stage of the outburst (e.g., rise, peak, decay, etc.; see lfT4l |9T1 T262I). 
There is a very strong linear anti-correlation between the fractional amplitude of the 
first overtone and the X-ray flux (something very similar was also observed in XTE 
J 1807-294; 11261 ). suggesting that the secondary hot spot becomes visible or more 
prominent as the accretion process becomes less intense. This may be explained if 
the inner disk radius moves in or out with respect to the accretion rate, rendering the 
secondary hot-spot visible lfT4l |9T i . If the inner disk radius r,„ is assumed to vary at 
different outburst stages, then the pulse profile changes of SAX J1808. 4-3658 can 
be explained in this way |263|. Variations in the disk-magnetosphere coupling may 
also explain explain some of the pulse shape and phase variations 112621 . Calculation 
of the light-curve of accreting NSs in global 3D-MHD simulations showed that if 
an octupolar magnetic field dominates the field configuration, double peaked pulse 
profiles can be reproduced H2411 . Similar pulse shape variations were seen to corre- 
late with the outburst stage of Swift J1756.9-2508 |26) IGR J1751 1-3057 (H and 
to some extent XTE J 1807-294 11321 . However, to date it is still unclear why pulse 
profiles change abruptly and unpredictably in some sources during certain outbursts 
but at other times have limited variability. 
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6 Long Term Evolution and Pulse Formation Process 



If AMXP pulsations are observed in different outbursts one can follow the long- 
term spin and orbital evolution. To date only five AMXPs have been monitored with 
high time resolution instruments in different outbursts: SAX J1808. 4-3658, IGR 
J00291+5934, XTE J1751-305, Swift J1756.9-2508 and NGC6440 X-2 (although 
with relatively low S/N and short outburst duration, it is difficult to constrain the 
long-term evolution of the latter two). 



6.1 Specific sources 



SAX J1808.4-3658 The best constrained AMXP is SAX J 1808.4-3658, for which 
secular spin evolution and orbital period variation have now been measured over 
a 13 year baseline lfl3l [T4l [T5l . As discussed previously, accretion torques have a 
small effect on the spin of SAX J1808. 4-3658 and have never been detected during 
any of the six outbursts monitored by RXTE. However, if one compares the constant 
spin frequency measured in each outburst, SAX J1808. 4-3658 is clearly spinning 
down at a constant rate (Bottom panel in Figure [8] [13]). The stability of the spin- 
down rate suggests that its origin is unrelated to propeller onset, since spin-down 
would then depend on the amount of mass propelled and on the duration of this 
phase, which can vary from outburst to outburst. Loss of angular momentum via 
emission of gravitational waves was also considered [14], and is interesting since 
this might also balance accretion torques, explaining the lack of measurable spin- 
up during the outbursts. However, fine tuning is also required in this case since the 
mechanism producing gravitational waves has to be triggered for the right amount 
of time so as to guarantee a constant spin-down between outbursts [13, 233 15 1. The 
most likely scenario appears to be a loss of angular momentum via magnetic-dipole 
radiation. This is expected for any rapidly rotating NS with a magnetic field, even if 
the NS does not turn on as a radio pulsar. The measured spin-down is of the order of 
- 1 1 5 Hz s "~ 1 and is consistent with a polar magnetic field of 2 x 1 8 G |T3] US El . 

The long-term orbital evolution of SAX J1808.4-3658 is also exciting. Theory 
states this should be driven by loss of angular momentum via emission of grav- 
itational waves from the tight binary orbit. For SAX J1808. 4-3658, the resulting 
orbital evolution should proceed on a timescale ~ 10 9 years (12641 : 

nm (^iVS+M 2 )' /3 D 8/3^ 

T GW =0.01 — P b ' Gyr (20) 

M NS M 2 b 

where is in hours and M^s and M 2 are the NS and donor masses. The observed 
value, however, implies an expansion of the orbit on a timescale of ~ 70 Myr |[T4l[T5l 
12651 [266 1. Further observations in 2011 provided a long enough baseline to detect 
an orbital period acceleration (Top panel in Figure [8] [13]). SAX J1808. 4-3658 
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Fig. 8 Orbital (top panel) 
and spin frequency evolu- 
tion (bottom panel) of SAX 
J1808.4-3658 over a baseline 
of 13 years. The top panel 
shows the delays accumulated 
by the pulsar as it passes 
through the ascending node. 
The data were consistent until 
2008 with a parabolic in- 
crease which is equivalent to 
a steady widening of the orbit 
(dotted line). The 2011 data 
revealed instead an accelera- 
tion of the orbit, as evidenced 
by the solid line in the top 
panel. The spin frequency has 
shown a steady decrease in 
between outbursts (bottom 
panel; Vo is a frequency offset 
of 400.975210 Hz), com- 
patible with a spin down of 
the order of -10~ 15 Hzs~'. 
(Figure from lfl3l ). 

may therefore behave like the "black widow" binary millisecond pulsars observed at 
radio frequencies 02671 with a short-term exchange of angular momentum between 
the donor star and the orbit. A similar suggestion was also previously proposed for 
this AMXP 1 266 265 1, invoking a "hidden" black widow scenario in which the 
AMXP switches on as a radio pulsar during quiescence and ablates its donor. 

IGR J00291+5934 Long-term spin down has also been measured between the 
2004 and 2008 outbursts of IGR J0029 1+5934. The AMXP rotation behaves like a 
saw-tooth, spinning up during the 2004 outburst then spinning down until the 2008 
outburst. The spin-down rate is sa — 3 x 10~ 15 Hzs~' |26, 33, 34], comparable to that 
observed in SAX J1808.4-3658. Although spin down between outbursts has been 
observed only once and it is not possible (until the next outburst) to verify whether it 
is constant, the most likely explanation also involves magnetic dipole radiation spin- 
down. If this saw-tooth behaviour is a good indicator of the secular spin evolution 
of this AMXP, then the net spin frequency is increasing at ss +2.5 x 10 _15 Hzs _1 . 
Indeed, the spin frequency at the beginning of the 2004 outburst was less than that 
observed in the 2008 outburst: the 2004 spin-up accelerated the pulsar more than it 
subsequently decelerated in quiescence. The relatively small baseline did not permit 
any useful constraints on orbital period evolution: this should be possible after the 
observation of the next outburst. 

XTE J1751-305 Comparing the spin frequencies from 2002, 2005, 2007 and 2009 
reveals a secular spin down of -5.5 x lO^Hzs- 1 |268|. This is again consistent 
with magnetic dipole spin-down. If the spin-up measured during the 2005 outburst 
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reported in |[T9l is not strongly affected by timing noise, then its behaviour resem- 
bles IGR J0029 1+5934, with spin frequency increasing in outburst and decreasing 
between outbursts in a saw-tooth. The inferred magnetic field is ~ 4 x 10 8 G, in 
line with the other two AMXPs constrained in this way. The low quality of the data 
taken in the 2002 outburst does not permit stringent constraints on orbital evolution. 

Swift J1756.9-2508 Observations during the outbursts in 2007 and 2009 give upper 
limits of |v| < 2 x 10~ 15 Hzs~' on the secular spin evolution of this AMXP [26|. 
This requires a magnetic field < 9 x 10 8 G. Magnetic field estimates for all four of 
these AMXPs are given in Table l6.ll 



Source Spin-Down in Quiescence Magnetic Field Secular Evolution Reference 

[Hzs-'] [10 8 G] 

SAX J 1808.4-3658 10 ls 1.5-2.5 Spin-down I13II14II151 

IGRJ00291+5934 -3 x lO" 15 1.5-2.0 Spin-up [26] 33 34 

XTEJ1751-305 -5.5xl0- 15 4 Spin-down [268] 

Swift J1756.9-2508 < |2| x 10~ 15 <9 ? [40] 

Table 4 Secular Evolution and Inferred Magnetic Fields in AMXPs. The magnetic field is 
inferred from the spin down observed in quiescence, and refers to field at the poles of the NS for a 
pure dipolar configuration and R = 10 km. 



6.2 The Maximum Spin Frequency of Neutron Stars 

None of the AMXPs discussed above has a spin rate that increases on long timescales, 
except for IGR J0029 1+5934. Even this source has a net acceleration so small that 
its spin frequency will change significantly only on timescales of several billion 
years. This raises the question on whether this behavior is the norm for AMXPs. 

For all realistic EoS of ultra-dense matter, NSs are stable at spin frequencies well 
in excess of 1000 Hz. The break-up frequency is well approximated by 

V max = 1230f^V / V^y 3/2 H Z (21) 

(where Mus and R^s refer to the non-rotating mass and radius of the NS under con- 
sideration, expression valid for arbitrary NS mass and EoS as long as the mass is not 
too close to the maximum permitted for that EoS [269, 270|). However the distri- 
bution of spin frequencies of the ensemble of AMXPs and NXPs (nuclear powered 
X-ray pulsars, see Introduction and Section 7) has an abrupt cutoff at about 730 Hz. 
This was first noticed in 2003 [92] and has been confirmed in later works with larger 
sample size B271I [32ll . Figure [9] shows the current distribution of spin frequencies 
for the 14 AMXPs and the 10 NXPs known (see Table 1 in [272 1 for a complete list 
of NXPs). So far no AMXP or radio millisecond pulsar has been found above this 
cutoff (the fastest radio millisecond pulsar has a spin of 716 Hz [273 1). So not only 
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Fig. 9 Histogram of the spin frequency of AMXPs and NXPs. The histogram is empty for fre- 
quencies larger than about 700 Hz. This is similar to what is found in radio pulsar data, where the 
sample size is considerably larger (with hundreds of millisecond radio pulsars). 



there are no pulsars with v 4 > 700 Hz, but also at least four AMXPs that should be 
accelerating in response to accretion are instead either decelerating or spinning up 
very slowly, on timescales of billions of years. 

Gravitational waves have been invoked to explain this cutoff, but as shown re- 
cently ||33l [34l I2331 , this cannot be the explanation for all AMXPs as it would 
require substantial fine-tuning to explain the X-ray timing observations. The cut- 
off might instead be related to the magnetic field evolution of NSs [223]. However, 
further investigation is needed to assess this fascinating question in a robust way 
since it is still unclear how the magnetic field of NSs evolves in response to accre- 
tion (274] [275] [276J \T77\ [278] [279). 



6.3 Why do most Low Mass X-ray Binaries not pulsate? 

Even after 15 years of high time resolution X-ray observations, only a few LMXBs 
have shown millisecond pulsations. Several LMXBs have pulsations at long periods 
(Table Q~|i but the number of pulsating systems is still small compared to the entire 
LMXB population. Several mechanisms have been proposed to solve this problem, 
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including burial of the magnetic field by accretion 1280, 12741 12811 , smearing of 
the pulsations by an optically thick corona [ 282 283, 284], smearing of pulsations 
due to gravitational light bending [285, 286 1, alignment of the NS magnetic and 
rotational axes 1 240 238 1 and onset of MHD instabilities at the disk/magnetospheric 
boundary [287 1. None of these models have yet been confirmed, although several 
new findings have helped to test their predictions. 

Magnetic Field Burial In this model, fresh unmagnetized accreted plasma dia- 
magnetically screens the NS magnetic field as it is deposited on the surface [280 
12741 [28 11. The field then re-emerges via Ohmic diffusion. Competition between the 
Ohmic diffusion and accretion timescale sets the behaviour of the magnetic field: if 
the accretion timescale is shorter than the Ohmic diffusion timescale, the magnetic 
field is buried deep in the ocean and outer crust. The two timescales can be defined 
as: 

^acc = ~r (22) 
m 

H 2 

Tdiff = — (23) 
*7 

where y = pdz is the column density, z is the direction orthogonal to the surface 
in the plane parallel geometry approximation, m is the local mass accretion rate 
per unit area (in units of gcm _2 s _1 ), H — ^r- is the pressure (defined as P) 
scale-height, and 77 is the magnetic diffusivity coefficient. This model predicts that 
LMXBs with very low mass accretion rates should be AMXPs, and that the most 
luminous LMXBs should not pulsate. This seems to be the case, since persistently 
pulsating AMXPs are indeed all faint LMXBs. However, the discovery of intermit- 
tent pulsations as short as ~ 100 seconds in Aql X-l poses severe problems for the 
model since it is unclear how a buried B field could emerge for this brief episode. 
Another interesting observational result comes from HETE J 1900. 1-2455, which, 
with a mean bolometric luminosity of 4.4 x 10 36 erg s -1 , has a mean mass accretion 
rate M rj 8 x lO _lo M yr _1 for an assumed accretion efficiency of 20%. This ac- 
cretion rate is about twice that of most AMXPs, but is smaller than the bolometric 
luminosity of SAX J1808.4-3658 during observations taken close to the peak of the 
2008 outburst [262|. This might indicate that the burial timescale is long compared 
to the outburst duration of most AMXPs. In HETE J1900. 1-2455 observational evi- 
dence for magnetic field burial has been found ll62l . with the magnetic field possibly 
reducing by almost three orders of magnitude on a timescale of about 100 days. If 
confirmed this indicates that indeed most LMXBs might be not pulsating because 
of the lack of an extended magnetosphere. Despite still open problems, the burial 
model remains therefore one of the most promising although the 1 -D nature of cur- 
rent models prevents the formulation of more sophisticated predictions and tests. 

Smearing from an Optically Thick Corona This model suggests that X-ray emis- 
sion from AMXPs is formed primarily via Comptonization processes in a relatively 
optically thin medium. If the Compton cloud is thick, then pulsations can be washed 
out. This scenario has been criticized on the grounds that spectral analyses of most 
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LMXBs show small optical depths, well below the limit required to smear out pul- 
sations 112881 , which should therefore be visible (but see [284]). It is not clear why 
in this model the intermittent AMXP SAX J1748.9-2021 shows pulsations that ap- 
pear and disappear on timescales of a few hundred seconds PT1 I40 I. uncorrected 
with spectral changes. However it provides a reasonable explanation for the time 
lags and pulse energy dependence, especially for the pulses of IGR J0029 1+5934 
which have a turnover at 7 keV that is otherwise difficult to explain [30, 145 1. 

Gravitational Light Bending For large stellar compactness (M/R), gravitational 
light bending can distort the path of X-ray photons up to the point where pulsations 
are strongly reduced [286. 289. 285 1. Persistent sources may therefore be less likely 
to show pulsations because they accrete more mass than transient systems. However, 
the total amount of mass that a NS has accreted does not depend on the current mass 
accretion rate and on the persistence of the X-ray flux, but on the evolutionary stage 
of the donor star and on the average mass transfer rate from the donor towards the 
accretor. This model also requires that the NSs in most LMXBs are more massive 
that those in AMXPs. Although possible, this is difficult to reconcile with the fact 
that most AMXPs have probably accreted for a long time (see for example [10|). 
The model also requires some special configuration of the accretion column, which 
existed for only a very short time and has never recurred, to explain the short pulsed 
episode of Aql X-l. Finally, the light bending effect never suppresses completely 
oscillation amplitudes since the Doppler shift due to the rapid NS rotation introduces 
anisotropies in the emission pattern an thus pulsations 12901 12531 . 

Magnetic and Rotational Axes alignment The possibility that the magnetic poles 
of NS migrate was first suggested by Ruderman in 1991 [240 1 . Vortex/fluxoid cou- 
pling in the NS interior links the spin and magnetic axes. Once the NS is spun-up by 
accretion, crustal shear stress breaks the crust into "plates" that drift towards the ro- 
tational poles. The phenomenology of an AMXP with a nearly aligned magnetic and 
rotational axis was developed recently [238] and explains fairly well the observed 
correlations between pulse amplitude, pulse phase and X-ray flux. The model can 
explain intermittency |291|, but implies that the persistent AMXPs should also 
show occasional "oscillation dropouts", a sudden and brief disappearance of pulsa- 
tions due to sporadic alignment of the hot spot and rotational axis. Such dropouts 
have not been observed, although the timescale over which this might happen may 
be too short compared to timescale over which the pulsations are measured (i.e., 
hundreds of seconds). The model also predicts correlations between excess back- 
ground noise produced by the hot spot motion which is anti-correlated with the 
pulse amplitude. This prediction could in principle be easily verified in the future 
with observational tests. 

MHD Instabilities Accretion onto a magnetized NS can proceed via channeled 
accretion but also via the Schwarzschild-Kruskal instability (the magnetic version 
of the Rayleigh-Taylor instability) in the disk equatorial plane ([292 293 1). This 
latter mode of accretion has now been observed in 3D-MHD simulations ( B294II ) 
and predicts a decoherence of the X-ray signal at high accretion rates, when several 
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"tongues" of plasma penetrate the magnetosphere and impact the NS surface at ran- 
dom positions. This model explains why bright LMXBs do not pulsate, and might 
explain why intermittent sources are on average brighter than persistent AMXPs. 
However, the intermittent AMXP SAX J1748. 9-2021 has shown pulsations that ap- 
pear and disappear in a broad range of accretion rates [HOI (rather than at a sharp 
threshold). It is also hard to explain why tongues never develop in any other AMXP 
when the outbursts reach peak luminosity, during which the accretion rates are com- 
parable or higher than observed in intermittent AMXPs. 



7 Thermonuclear Bursts 



Source 


Burst oscillations Phase locking 


References 


Accreting Millisecond X-Ray Pulsars 


SAX J1808. 4-3658 


Yes 


No 


74. 75 92 295 82 


29611297 


XTE J1814— 338 


Yes 


Yes 


I298||299||300||301 


236 297 


Aql X-l 


Yes 


No 


I302ll303ll37)51ll60ll305l 


3rj6ll2?7 


SAX J1748.9-2021 


No 




169 307 41 297 


HETE J 1900. 1-2455 


Yes 


No 


308 144 61 


297 309 


IGR J 175 11 -3057 


Yes 


No 


[182 183 310 44 


IGR J 17498-2921 


Yes 


Unknown 




20311204 


Swift J 1749.4-2807 


No 






19lll200 


Mildly Recycled X-ray Pulsar 


IGR J 17480-2466 


Yes 


Yes 


49 311 312 313 314 315 



Table 5 Burst oscillation properties for the eight AMXPs that have shown thermonuclear 
bursts. The table indicates whether a source shows burst oscillations and whether their phases 
match and follow the evolution of the accretion-powered pulses observed before and after the 
bursts ("Phase Locking"). 



Type I X-ray bursts are thermonuclear explosions triggered by unstable burn- 
ing of hydrogen or helium on the layers that build up on the surfaces of accret- 
ing NSs. The basic cause is an imbalance between nuclear heating and radiative 
cooling, which leads to the runaway of temperature-dependent thermonuclear reac- 
tions 1 304]. There are however many aspects of burst physics that remain puzzling, 
such as discrepancies between predicted and observed recurrence times, and the 
phenomenon of burst oscillations 03161 1272 1. The latter are high frequency oscil- 
lations (11-600 Hz) often detected in the X-ray flux produced by thermonuclear 
bursts 12971 , and whose cause remains unknown. AMXPs can offer particular in- 
sight into burst physics and the burst oscillation mechanism, for two reasons. Firstly, 
they are the only sources where we have an independent measurement of the stel- 
lar spin, allowing us to measure the dependence on rotation. Secondly, they are the 
only sources where we can make a reasonable assessment of the effects of the mag- 
netic field, and the role of uneven fuel distribution (due to magnetic channeling of 
material onto the poles of the star). 
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Table 5 summarizes the AMXPs that have exhibited bursts, and which of these 
have shown burst oscillations. Bursts are expected and observed to occur only within 
a range of accretion rates from ~ 2 — 30% of the Eddington rate 03171 13181 13191 , 
and in this regard the occurrence of bursts in the AMXP population accords with 
the behaviour seen in the non-pulsing bursters. The sources that do not show bursts 
are the five ultracompact sources (which are thought to have rather low accretion 
rates [316|), likely below the minimum required for regular bursting [320]) and 
IGR J0029 1+5934, which is also thought to have a rather low accretion rate [29 1 . 
None of the AMXPs has yet shown a super-burst, longer duration bursts thought to 
be due to unstable carbon burning 113211 . 

Burst properties. Bursting behaviour is known to be extremely variable, and the 
bursts of the AMXPs are no exception to this rule. To date there has not been a 
comprehensive study assessing whether the bursts from the AMXPs have properties 
that differ in any statistically significant way from the non-pulsing sources. How- 
ever just as for the non-pulsing bursters, bursts from the AMXPs show variable 
recurrence times, peak fluxes and fiuences, and durations, even when accretion rate 
seems to be relatively stable. There are nonetheless some features of AMXP bursts 
that are worthy of note, in particular when they are compared with bursting sources 
with higher magnetic fields like the mildly recycled pulsar IGR J17480-2466. 

Magnetic fields and bursting. Bursts from AMXPs confirm that magnetic fields 
~10 8 -10 10 G(Table[D are no impediment to thermonuclear bursting. 

Burst recurrence times. Many sources show burst recurrence times of less than an 
hour, something that is problematic for theoretical models 13221 . The only AMXP 
in this short recurrence time group is the highly intermittent source Aql X-l. By 
contrast the mildly recycled pulsar IGR J 17480-2466 has the shortest recurrence 
time between thermonuclear bursts yet recorded, at 3.3 minutes A3 1 II . a possible 
sign of confined burning. 

Rotation and bursting. The AMXPs, which have independent estimates of spin rate, 
confirm that bursting can occur for spin rates of up to a few hundred Hz (although the 
fastest rotating burster, 4U 1608-522, is a burst oscillation source not an accretion- 
powered pulsar, so the spin measure is more indirect [323 1). By contrast the mildly- 
recycled accreting pulsar IGR J 17480-2466 is the bursting source with burst oscil- 
lations with the lowest measured rotation rate, at 1 1 Hz. For a star rotating at ~ 10 
Hz, the Coriolis force is not dynamically relevant 1149] 12721 . indicating that Coriolis 
force induced confinement of the igniting patch is not essential to the development 
of X-ray bursts [ 324 1 or to very short burst recurrence times (see previous point). 

Cooling as a signature of thermonuclear bursts. Cooling during the tails of bursts 
(inferred from blackbody fits to the burst spectra) has long been regarded as a sig- 
nature trait of thermonuclear Type I (rather than accretion -powered Type II) bursts. 
Detailed studies of the bursts from the mildly recycled pulsar IGR J 1 7480-2466 have 
now shown that cooling is a sufficient but not necessary condition for the identifica- 
tion of thermonuclear bursts: under certain circumstances, particularly at high accre- 
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tion rates, the cooling signature may not be readily detectable 113 1 21 13 1 3l 13 141 |315l . 
Bursts from the AMXPs do however show cooling. 

Marginally stable burning. The mildly recycled pulsar IGR J 17480-2466 shows 
mHz QPOs that are interpreted as marginally stable burning of H and He on the NS 
surface [325 3151 . The mHz QPOs are observed when the accretion rate rises and 
the source transitions from unstable burning (X-ray bursts) to stable nuclear burning 
(no bursts). The only AMXP to have shown mHz QPOs is the highest accretion rate 
source, the intermittent pulsar Aql X-l 113261 . 

Pure helium bursts. SAX J 1808.4-3658 has very bright, short bursts that are thought 
to be the only secure example of bursting occurring in a layer of pure helium, for 
a source where the accreted fuel is a mix of hydrogen and helium [82|. The only 
other sources that are thought to have pure helium bursts are ultra-compact sources 
where the donor star is too small for there to be significant hydrogen in the accreted 
material. Pure helium bursting is thought to occur in mixed fuel sources only in a 
very narrow range of accretion rates where hydrogen can burn stably to form helium 
before the helium ignites unstably 03171 . 

Bursting and intermittent accretion-powered pulsations. Two of the intermittent 
AMXPs, HETE J1900.1-2455 and SAX J1748.9-2021, show intermittent pulsations 
with variable amplitudes at times when the source is also bursting IBTl |4D . There 
does not appear to be any clear-cut causal relationship between the two phenomena 
(and the third intermittent AMXP, Aql X-l, does not show bursts at the time of its 
one intermittent pulsation episode, [38 1). Nonetheless there has been speculation in 
the literature as to whether the two might be related. Radiation from the burst might 
perhaps lead to a temporary disruption of the inner edge of the disk, perturbing the 
flow of material onto the magnetic poles and moving the accretion footprint to a 
location more favorable for the formation of pulsations H291I . 

Burst oscillations Where the AMXPs have had the biggest impact is in our under- 
standing of burst oscillations (see [272 1 for a more comprehensive review of this 
topic). Burst oscillations were first observed in 1996, from the non-pulsar 4U 1728- 
34 03271 . The oscillations manifested as a coherent signal at about 363 Hz, with an 
upwards drift in frequency of ~ 1 Hz as the bursts progressed, towards an asymptotic 
maximum in the tails. The fact that the same frequency was seen in multiple bursts 
from the same source, and the stability of the asymptotic maximum frequency, sug- 
gested a link to a stable clock such as the stellar spin rate [328 1. Firm identification 
of burst oscillation frequency with the spin frequency however had to wait until the 
first robust detection of burst oscillations from an AMXP [92 1. Although the two 
frequencies appear to be very close (within a few Hz), however, the separation does 
differ markedly for different AMXPs. The AMXP XTE J1814-338 and the mildly 
recycled pulsar IGR J 17480-2466 have burst oscillation frequencies that agree with 
the spin frequency very closely (within the error bars, - 10~ 8 Hz for XTE J1814- 
338 j298H3UQl|236l, - 10~ 4 Hz for IGR J17480-2466 |49)). SAX J1808.4-3658 
and IGR J 175 11-3057, by contrast, have burst oscillations that show a rapid in- 
crease in frequency in the burst rise, slightly overshooting the spin frequency, and 
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then stabilizing to within ~ 10 mHz of the spin frequency in the tails |92l 11821 . 
The intermittent pulsars Aql X-l and HETE J 1900. 1-2455 show burst oscillations 
with slow drifts throughout the bursts, more similar to those of the non-pulsars, and 
in these cases the asymptotic maximum frequency is ~ 1 Hz below the spin fre- 
quency [160, 329 38, 309 1. Any model for the burst oscillation mechanism must 
be able to explain this diversity in frequency drift and the small variations in offset 
from the spin frequency. 

There are interesting differences in the properties of the burst oscillations from 
the pulsars as compared to those of the non-pulsars. For the non-pulsing bursters, 
burst oscillations tend only to be seen when the source is in the soft (high accretion 
rate) state [297 1. The two intermittent pulsars conform to this rule: indeed HETE 
J 1900. 1-2455 has shown burst oscillations in only one burst, when the source en- 
tered an unusually soft state [ 309] . By contrast four of the five persistent pulsars with 
bursts have had burst oscillations in all of their bursts, even though these sources 
tend to be in the hard (low accretion rate) state. For the fifth, IGR J17498-2921, the 
data are of much poorer quality: oscillations are seen only in the brightest bursts but 
the upper limits on the presence of detections in the weaker bursts are comparable to 
the amplitudes detected in the brighter bursts [204] . The duration of burst oscillation 
trains also appears to differ. For SAX J1808.4-3658, XTE 1814-338 and the mildly 
recycled pulsar IGR J17480-2466, burst oscillations persist throughout the bursts 
(except during episodes of Photospheric Radius Expansion in the bright bursts from 
SAX 1808.4-3658). For the non-pulsars, although burst oscillations are sometimes 
detected throughout bursts they are more commonly detected in the tails [297 1, and 
this is also the case for the intermittent accretion-powered pulsars. In this regard the 
persistent pulsar IGR J 1751 1-3057 seems to be a transitional object: at lower accre- 
tion rates oscillations are detectable throughout bursts, but as accretion rate rises the 
burst oscillation signal vanishes from the burst rise [ 1 82 1 . Burst oscillations from 
the pulsars have higher harmonic content than burst oscillations from the intermit- 
tent pulsars and non-pulsars [123T1 [92l l298l [3001 1309J II 82l l49ll . They also have rather 
different amplitude-energy relations. Burst oscillations from the persistent pulsars 
show the same amplitude-energy relationship as their accretion-powered pulsations, 
irrespective of whether the latter rise or fall with energy but with the latter being 
more common (fl92l|298l[309l[l82l|49l|48l|^ and see also SectionED. The inter- 
mittent pulsars and non-pulsars, by contrast, have burst oscillation amplitudes that 
rise with energy 113301 1 1 84l . The cause of these various differences in burst oscilla- 
tion properties between the pulsars and non-pulsars remains as yet unclear, but the 
most obvious hypothesis is that the differences are due to the effects of a dynami- 
cally important magnetic field. 

One can also use the pulsars to compare the properties of the accretion-powered 
pulsations (where the peak of the emission is presumed to be centered on the mag- 
netic poles) and their burst oscillations. Burst oscillation amplitudes vary substan- 
tially from source to source, but are in all cases within a few percent of (and most 
lower than) the amplitude of the accretion-powered pulsations at the time of the 
burst []92l [3001 [3091 [T82l |49l HH . Pulsar burst oscillations show no statistically 
significant phase lags as a function of energy If330l [30T| [309), unlike accretion- 
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powered pulsations which show soft lags (1260, [33T1 |2T1 332 333, 1 8 1 1 ; see also 
Section IBTTT i. Of particular interest is the strong phase-locking between accretion- 
powered pulsations and burst oscillations seen in XTE J18 14-338 and IGR J 17480- 
2466 [236, 49 1, with very small phase offset ( < 3° and < 10° respectively), imply- 
ing that the hot spots responsible for the two sets of pulsations are longitudinally 
coincident in these two sources. 

So what do these results tell us about the burst oscillation mechanism ? Although 
the mechanism is not yet understood, substantial progress has been made and the 
accretion-powered pulsars have played a key role. Their primary contribution, of 
course, has been to highlight the very close relationship between the burst oscilla- 
tion frequency and the spin frequency of the star. This fact underpins the two main 
classes of model: hot spot models and global mode models. Small temperature vari- 
ations in the surface layers mean that ignition is expected to begin at a point, with 
a flame front then spreading out across the star [334 1 . This flame may then either 
stall, confining the burning to a small region (the hot spot model) or excite large- 
scale waves in the surface layers (global mode models). The resulting temperature 
asymmetry gives rise to the burst oscillations. In hot spot models, the relationship 
with the spin rate is straightforward since the hot spot should be near stationary 
in the rotating frame of the star 113271 . The main open theoretical question for hot 
spot models is what might cause the flame to stall, with both magnetic and Coriolis 
forces under consideration [335, 324, 49 1. For global mode models, the requirement 
that burst oscillation frequency be very close to spin frequency puts very stringent 
restrictions on the types of surface modes that might be responsible [336, 337, 338 1. 
Efforts are now ongoing to develop both classes of model to determine whether they 
could match the rest of the observed properties [272]. The differences between the 
burst oscillations from the pulsars and the non-pulsars are informative, and suggest 
a role for the magnetic field. What is still unclear, however, is whether we are see- 
ing one mechanism with a continuum of properties set by stellar parameters such 
as the magnetic field - or whether two different mechanisms may be required. The 
apparent gradual changes in burst oscillation properties from the non-pulsars and 
intermittent pulsars (where magnetic field is likely to be weak) to the persistent 
pulsars (with stronger fields) had until recently favoured the former. However new 
results from the pulsars point increasingly towards the possibility that we are in- 
deed seeing two different mechanisms. The development of burst oscillations in the 
slowly rotating pulsar IGR J 17480-2466 is particularly hard to understand in the 
context of global modes or Coriolis force induced hot spots [ 49 1 . Magnetic confine- 
ment of the flame front is at present the only plausible explanation for this source, 
and would also explain the extraordinary phase-locking between accretion-powered 
pulsations and burst oscillations. This mechanism might also operate in the other 
strong magnetic field source with phase-locked burst oscillations, the AMXP XTE 
J18 14-338 1 236, 49 1. However for the other sources, with weaker magnetic fields, 
magnetic confinement is unlikely to be effective. For these sources global mode 
models remain a good possibility, with the differences in the pulsar burst oscilla- 
tions being perhaps due to magnetic modifications to the mode structure. 
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8 Aperiodic Variability and kHz QPOs 

The study of aperiodic signals in power spectra of AMXPs reveals a rich phe- 
nomenology which has helped us to understand the physics of accretion disks and 
compact objects. All AMXPs are "atoll sources", so-called because of their pattern 
in the X-ray colour-colour diagram. Atoll sources are relatively low/intermediate- 
luminosity LMXBs and are believed to possess low magnetic fields compared to the 
brightest LMXBs (the "Z sources"). The aperiodic variability of AMXPs includes 
phenomena at very high frequencies (the kHz QPOs) but also at lower frequency, 
such as the strong 1 Hz modulation observed in SAX J1808.4-3658 and NGC6440 
X-2 (Sections 13.11 and [3 . 1 1 \ . Here we highlight the main aperiodic phenomena pe- 
culiar to AMXPs and refer to Chapter 1 of this book for a more detailed discussion. 

The fastest variability in LMXBs is observed as kHz QPOs, which appear at 
frequencies of up to ~ 1300 Hz. Sometimes two kHz QPOs are observed simulta- 
neously as "twins", with an "upper" QPO at frequency v„ and a "lower" QPO at 
frequency V/. Such fast variability is observed only in LMXBs with NS accretors, 
and is considered to be a signature of the presence of a NS in the absence of pulsa- 
tions or thermonuclear bursts. All models that try to explain kHz QPOs associate at 
least one of the frequencies with the orbital motion of plasma in the accretion disk. 
Such high frequency variability is expected given the short dynamical timescale 
(0.1-1 ms) characterizing the motion of matter near NSs. Such rapid variability is 
never observed in black hole LMXBsH possibly because they have larger mass and 
hence a larger innermost stable circular orbit (ISCO). The ISCO is a limiting radius 
(at about 6GM/c 2 for non-rotating NSs) set by General Relativity that defines the 
the smallest stable orbit in an accretion disk. The frequency of kHz QPOs is close 
to the value expected for stable Keplerian orbits around a NS, and can therefore 
be used to place constraints on the mass and radius of the NS (and hence the EoS 
of ultra-dense matter). The size of the orbit where kHz QPOs are produced must be 
larger than both NS radius and ISCO, both of which depend on NS mass (FigurefTOii. 

The four AMXPs with twin kHz QPOs (SAX J1808.4-3658, XTE J1807-294, 
IGR J 1 75 1 1-305 and Aql X- 1 ; Section© have been used to test kHz QPO formation 
models (e.g., the rejection of the beat-frequency model; Section [3~TT i and have illu- 
minated the relation with the NS spin. Taking into account also the NXPs with twin 
kHz QPOs, it is clear that the relation between kHz QPO separation 4v = V„-V/ 
and spin frequency V s is not as simple as originally thought: Av/v s is not always 
equal to 1 or 0.5. Av and may even be unrelated: instead we see an average 
(Av) « 300 Hz for all sources except AMXPs J340l[34j]- Only for the AMXPs, 
where the magnetic field is strong enough to channel accretion, is the ratio Av/v s 
either 1 or 0.5. Some predictions of the sonic -point spin resonance and relativistic- 
resonance models B3421 1343 1 are consistent with the AMXP results. In sonic-point 
spin resonance models, the magnetic and radiation fields rotating with the star ex- 
cite a vertical motion of the disk gas at the "spin-resonance" radius r sr . There, the 
vertical epicyclic frequency equals the difference between the orbital frequency and 



3 Fast variability appears in black hole LMXBs below ~ 500 Hz as "high frequency QPOs". 
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Fig. 10 Constraints on the 
Mass-Radius relation of NSs 
from kHz QPOs. The solid 
curves refer to a sample of 
different EoS models (for 
non-rotating NSs), and the 
thick solid and dashed lines 
enclose regions of the M-R 
diagram compatible with the 
observation of kHz QPOs at 
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the spin. Depending on the "dumpiness" of the flow at r rs , the X-ray flux exhibits 
A v equal to either V s (smooth flow) or 0.5 Vj (clumpy flow). In relativistic resonance 
models, a non-linear 1:2 or 1:3 resonance between orbital and radial epicyclic mo- 
tion emerges as a consequence of the deviation of the strong gravitational potential 
from the Newtonian 1 jr. A third model suggests that kHz QPOs are related to the 
three general relativistic epicyclic frequencies 113441 or to the precession of frame 
dragging (Lense-Thirring precession 113451 ). So far none of these models satisfacto- 
rily accounts for all of the rich phenomenology of kHz QPOs. 

Determining the origin of QPOs, in particular kHz QPOs is of fundamental im- 
portance as this is a phenomenon that is common in LMXBs. If a clear relation 
with the NS mass and/or radius is confirmed, it has the potential to place severe 
constraints on the EoS of ultra-dense matter. This field has however seen few ad- 
vancements in the past few years, at least on the theoretical side. New attempts are 
now using numerical simulations to interpret the observations and test some of the 
original models [346 347, 348, 237 1 . Most of these simulations use black hole ac- 
cretors (but see [237 1) which do not have kHz QPOs, but the results are promising 
in terms of understanding the generation of variability in LMXBs as a whole. 



9 Open Problems and Final Remarks 

Since the first discovery in 1998, the AMXPs have provided an incredibly body of 
observational data allowing us to understand NSs and their evolution. They have 
allowed the study of extreme phenomena on the surface of the NSs, in the strong 
gravity regime at relativistic rotational velocities, and have helped us to understand 
how radio pulsars form and evolve within the framework of the recycling scenario. 

Proposed X-ray missions like the Indian ASTROS AT, the ESA mission LOFT and 
the NASA AXTAR satellite will guarantee huge advances in our knowledge of these 
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systems, building on the extraordinary results from RXTE, XMM-Newton, Chandra 
and Swift. Long-term monitoring of future AMXP outbursts will allow measurement 
of long-term spin and orbital behaviour, something that has at present only been 
done for four systems. Increased collecting area will increase sensitivity to pulsa- 
tions by orders of magnitude, leading (we hope) to the discovery of sub-millisecond 
pulsars or strengthening the 730 Hz spin distribution cutoff. In the former case this 
would place very tight constraints on the EoS of ultra-dense matter. Modeling of 
high signal to noise pulse profiles will reveal the subtle features imparted by strong 
gravity, enabling us to measure the NS compactness. We will be able to test models 
for kHz QPOs and burst oscillations with sufficient precision that we can finally 
determine the underpinning mechanism. 

Of the many open problems, there are several that stand out. One big mystery is 
why do so very few LMXBs show pulsations? If magnetic field evolution and decay 
is involved, we may finally be able to understand the origin and the evolution of 
nature's strongest magnetic fields. What happens when accretion takes place during 
an outburst? Why we do not detect the effect of accretion torques in some AMXPs, 
contrary to expectations ? Are gravitational waves involved as a braking mechanism 
that counteracts the effect of accretion or are we again missing some fundamental 
ingredient? When do accreting pulsar switch on as radio pulsars ? We look forward 
to discovering the answers to these questions in the years ahead. 
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